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' We establish an explicit correspondence between perturbative and nonperturbative results in the 

problem of quantum decoherence in disordered conductors. We demonstrate that the dephasing 
time r v cannot be unambiguously extracted from a perturbative calculation. We show that the 
effect of the electron-electron interaction on the magnetoconductance is described by the function 
, Ad(t) exp( — fd(t)). The dephasing time is determined by fd(t), i.e. in order to evaluate t v it 

is sufficient to perform a nonperturbative analysis with an exponential accuracy. The effect of 
' interaction on the pre-exponent Ad(t) is important if one calculates the interaction-dependent part 

of the weak localization correction for strong magnetic fields. The zero temperature dephasing 
time drops out of this correction in the first order due to the exact cancellation of the linear in 
time T-independent contributions from the exponent and the pre-exponent. Nonlinear in time T- 
independent contributions do not cancel out already in the first order of the perturbation theory 
■ and yield an additional contribution to dephasing at all temperatures including T = 0. 

Jo 
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I. INTRODUCTION 

' Recent experiments by Mohanty, Jariwala and Webb§ strongly indicate an intrinsic nature of a low temperature 
saturation of the electron decoherence time t v in disordered conductorsou. It was argued in Refs. [l],^| that zero point 
fluctuations of electrons could be responsible for a finite dephasing at low temperatures. These as well as various other 
O , recent experimental results attract a lot of attention to the fundamental role of interactions in disordered mesoscopic 
systems. „ 
I A theory of the above phenomenon!!! was proposed in our papersOtl We demonstrated that electron-electron 
. interactions-ill disordered systems can indeed be responsible for a nonzero electron decoherence rate down to T = 0. 
Our resultaj'B are in a good agreement with experimental findings!!]. We also arguedQ that this interaction-induced 
decoherence, has the same physical nature as in the case of a quantum particle interacting with a bath of harmonic 
oscillatorstIM 

The low temperature saturation of the decoherence rate on a level predicted in Refs. has serious thmretical 
consequences. Therefore it is not surprising that .these predictions initiated intensive theoretical debatesEjllj. In 
Q\ \ contrast to Refs. [?], various authorsaO'q'EaOljj argued that interaction-induced electron dephasing at T — > 
is not possible. Vavilov and Ambegaokar (VA)ll3 argued the quantum correction to the classical resulfa'cl should be 
small at least in the limit Tt^ ^> 1. 

It should be emphasized that the above discussion goes far beyond the problem of electron dephasing only. This 
discussion is important for a general understanding of the role of the electron-electron interactions in mesoscopic 
systems at low temperatures. According e.g. to Aleiner, Altshuler and Gershenzon (AAG)I3 this role is merely to 
provide a (temperature dependent) renormalization of a disordered potential of impurities. Within this picture, at 
sufficiently low T (when the effect of thermal fluctuations is small and can be neglected) electrons propagate in an 
effective inhomogeneous static potential which should be determined self-consistently in the presence of Coulomb 
. ' interaction. If so, electron scattering on such astatic potential is not any different from that on static impurities and, 
hence, cannot lead to dephasing. Our resultsOu suggest a different picture, according to which dynamical effects are 
important at all temperatures down to T = and the high frequency "quantum" modes with w > T do contribute to 
dephasing. ^^^^^ 

Two main arguments supporting the first ("static") picture are usually discussedljOEjOllijlCh^ first argument 
is quite general and is not necessarily related to electrons in a disordered metal. One can argueE3EjllZl that a particle 
with energy ~ T cannot excite harmonic oscillators with frequencies lu > T and, hence the latter will at most lead to 
renormalization effects. It is easy to observe, however, that this argument explicitly contradicts to the exact results 
obtained e.g. within the Caldeira-Leggett modelO. There even at T = the off-diagonal elements of the particle 
density matrix decay at a finite length set by interaction. This effect is due to all high frequency modes of the 
effective environment, i.e. the picture is by no means "static" (see also Refs. 0,f2|JlJ for further discussion). 
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One can also modify the above argument and conjecturet3 that the system of electrons can behave differently from 
a bosonic oneu because of the Pauli principle which restricts scattering space for electrons at low T and, hence, their 
ability to exchange energy—Again, this argument contradicts to the well known results obtained for fermionic systems. 
E.g. it is well establishedEj'EIll that tunneling electrons exchange energy with the effective environment (formed by 
other electrons in the leads) even at T = 0. This exchange xesults in the temperature independent broadening of 
the effective-energy distribution P{E) for tunneling electronstj. This sOpGalled "P(E)-theoty" was verified in many 
experimental. A close formal and physical similarity between the theorylj and our analysisa is discussed in Ref. |2^. 

The second argument against the possibility of the interaction-induced saturation of t v is purely formal. It is based 
on a perturbative calculation by AAQI3. These authors claimedtllEj that the results of this calculation explicitly 
contradict to our resultsB'lj and, hence, the latter are inc orre ct. However, a convincing comparison between the two 
calculations was not presented. Furthermore, it our ReplyHLJ we pointed out that the origin of the above controversy 
lies deeper, and the AAG's suggestion that our calculation is "profoundly incorrect" can not be taken seriously. 
We arguecLthat both approaches do agree on a perturbative level and the key difference between them is that our 
calculations is nonperturbative while the analysis by AAG does not go beyond the first order in the interaction 
and, on top of that, involves additional approximatipas not contained in our papem. For instance, for the exactly 
solvable Caldeira-Leggett model we demonstratedlI3ii-J that within the perturbative approach involving analogous 
approximations one arrives at incorrect results and misses the effect of quantum decoherence at low temperatures. 

Motivated by this discussion as well as by the fundamental importance of the problem we have undertaken an 
additional analysis of the effect of interaction-induced deaoherence in disordered metals. This analysis will help us 
to demonstrate the actual relation between our approach^ and that of AAGlla. Since it is hardly possible to settle 
a calculational dispute without presenting sufficiently many details, in this paper we made an effort to provide the 
reader with such details of our calculation. 

The structure of the paper is as follows. In section 2 we will demonstrate a principal insufficiency of the perturbation 
theory in the interaction for the problem of quantum dephasing. We will argue that t v cannot be unambiguously 
extracted even from a correct perturbative calculation. In section 3 we extend our nonperturbative calculation!! We 
will carry out a complete analysis of the problem with the exponential accuracy. We will also present semi-quantitative 
arguments which, however, will be sufficient in order to understand the effect of interaction on the pre-exponent. In 
section 4 we perform a detailed perturbative calculation and demonstrate that at low T some previous perturbative 
results are based on several insufficient approximations, the main of which-is the golden rule approximation. We 
also establish an explicit relation between nonperturbativeQ and perturbativdlj calculations. A close formal similarity 
between the problem in question and the exactly solvable Caldeira-Leggett model is analyzed in section 5. In section 
6 we briefly summarize our main observations. For the sake of convenience we will briefly announce the main steps 
of our calculation in the beginning of each section Some further technical details are presented in Appendices A, B, 
C and D. In Appendix E we discuss the resultSujiiZl. 



II. INSUFFICIENCY OF THE PERTURBATION THEORY 



In this section we will demonstrate a principal insufficiency of a perturbative (in the interaction) approach to the 
problem of quantum dephasing. In the subsection A we will present some general remarks concerning the role of 
the perturbation theory for the problem of a quantum mechanical particle interacting with other quantum degrees of 
freedom. In the subsection B we discuss the relation between perturbative and nonperturbative calculations of the 
magnetoconductance and the decoherence time t v in disordered conductors. 



A. General remarks 



The time evolution of the density matrix of such a particle is defined by the following equation: 

p(t,X-Lf,X2f) = J dx u dx 2i J(t,X l f,X2f,Xli,X2i)p(0,Xu,X 2 i), (1) 

where x is the particle coordinate. The kernel J depends on the Feynman- Vernon influence functionallBEfl and contains 
the full information about the effect of interaction. This kernel can formally be expanded in powers of the interaction 
strength 

oo 

m = Y J j {n \t). (2) 

n=0 
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The "noninteracting" kernel does not change the state of the system (provided its initial state is an eigenstate 
of the noninteracting Hamiltonian) and in this sense it is equivalent to the unity operator. All other terms of this 
expansion grow with time the faster the larger the number n is. As a result in general all these terms (||) become 
important for sufficiently long times. Hence, the perturbation theory in the interaction (which amounts to keeping 
only several first terms of the expansion (^|)) is equivalent to the short time expansion of the exact density matrix. 
Thus in general this perturbation theory cannot correctly describe the long time behavior of the interacting system 
no matter how weak the interaction is. 

Keeping only the term in (||) under some additional assumptions one can express the probability Wa for the 
particle x to remain in its initial state as 

Wu(t) = 1 - f dh f 1 dt 2 K{t u t 2 ), (3) 
Jo Jo 

where the kernel K(t\,t2) can be derived from the influence functional and will not be specified here. In equilibrium 
one usually has K(t\, t^) = K{t\ — t<i). Eq. (J3J) applies at short times, when the second term is still much smaller than 
unity. But even in this limit the correct information can be missed by insufficient approximations. For instance, the 
frequently used approximation amounts to retaining only the u = term in the kernel K u — J (dr /2tt)K(t) exp(iwr). 
Within this so-called golden rule approximation one finds 

W u (t) = l-Tt, T = ttK lj=0 . (4) 

Furthermore, assuming that the effect of higher order terms in the expansion (|^) can be accounted for by exponentiating 
the last term in (^) one immediately arrives at Wu{t) = exp(— Ft). 

Obviously the above set of approximations is justified only in special cases. For instance, the golden rule approxima- 
tion can work only provided the kernel K(t) decays rapidly as compared to other relevant time scales in the problem. 
This could be the case e.g. at sufficiently high temperatures. In general, and especially in the low temperature limit, 
the golden rule approximation (^) fails. And it is particularly dangerous to combine the short time perturbative 
expansion with the long time golden rule approximation. E.g. if K u —o happens to be zero, it would follow from (^) 
that the particle will stay in its initial state forever even in the presence of interaction. Obviously this cannot be the 
case. The exponential decay of the probability Wu in time is also an artifact of the golden rule approximation. In 
general the time dynamics of an interacting system is much more complicated, and it should be determined from eq. 

©• n . 

In eq. (pi) is usually implied that the initial density matrix p(0) does not coincide with the exact reduced equilih ri j um 
density matrix for the interacting system. The standard approach is simply to factorize the initial density matrixEacJ, 
i.e. to represent it as a product of the particle density matrix p(0, xu, x 2 i) and the equilibrium density matrix of 
all other degrees of freedom. In this case, even if initially both the particle and the environment were in their 
noninteracting ground states at T = 0, the relaxation process occurs because the factorized density matrix does not 
describe the ground state of the interacting system. One could question the relevance of such initial conditions e.g. 
to the problem of electron transport in disordered conductors in the presence of interaction. Indeed, in this case the 
density matrix is never factorized and no time evolution can be expected for the equilibrium density matrix of the 
whole interacting system. Hence, at T = in equilibrium no relaxation should occur. 

In order to clarify the situation let recall the formal expression for the conductivity (see the eq. (|A15|)) 



croc J dt 1 J drudr 2i J(t - t' ,dr 1 f,dr 2 f,dr li ,dr 2i )p cff (r li ,r 2i ), (5) 

— oo 

where p c g (ru, r 2 i) = {ru — r 2 i)pa{ru, r 2 i) and po is the equilibrium electron density matrix. We observe that the 
effective initial density matrix p e g in this expression is strongly perturbed at all T as compared to po due to the factor 
(ru — r 2 i). Therefore relaxation always takes place in our problem. Since for a dissipative system relaxation times do 
not depend on the initial conditions, one can safely assume the initial density matrix to be, for instance, facto rized . 
Actually the same assumption is used within the diagrammatic approach: a complete equivalence between eq. ( |A15| ) 
(factorized density matrix) and the diagrammatic expression for the conductancet3 was demonstrated in Appendix A 
in the first order in the interaction. 



B. Magnetoconductance 



The weak localization correction to the conductivity of a disordered metal can be expressed in the following form 
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5a d {H) = J dtA d (t) exp (-t/r H - f d (t)) , (6) 

Te 

r e = 1/vf is the elastic electron mean free time. The function fd(t) increases with time and describes the Cooperon 
decay due to interaction (fd{t) equals to zero without interaction). The presence of the magnetic field H causes an 
additional decay on a time scale ~ th- By varying the magnetic field and thus th (which decreases with increasing 
H) one can extract information about the interaction-induced decoherence directly from the magnetoconductance 
measurements . 

The pre-exponential function Ad(t) without interaction is A d °\t) — l/(47rZ?i) d / 2 . In the presence of interaction 
the function Adit) will, of course, depend on the interaction as well. As it is demonstrated below, this dependence 
is can be ignored while calculating the decoherence time t v which should only be extracted from the function fd(t) 
in the exponent of @. This is the procedure of Ref. |[ However, the dependence of the pre-exponent Ad(t) on the 
interaction is important if one wants to recover the sublcading in th/t v term in the expression for Sad(H) in the 
limit of a strong magnetic field tjj <C t v . In this case only short times t < th contribute to the integral (|6|) and it is 
sufficient to perform a short time expansion of both exp(— f d (t)) and Ad(t). This expansion mixes terms important 
and unimportant for dephasing and in general makes it impossible to extract correct information about the dephasing 
time Ttp from the perturbation theory even in the limit of strong magnetic fields th <r r 

In order to illustrate this conclusion let us restrict ourselves to a quasi- Id case. The expression @ may then be 
rewritten as follows 

= "^F / ^ t/TH F{t/r v ), (7) 
o 

where the function F(t/r v ) accounts for the interaction. Note, that the function F can (and in general does) depend 
not only on one but on several parameters F = F(t/ri,t/T2, ...,t/r n ). In this section we will assume that F depends 
on only one parameter t v . This is sufficient for our purposes. 

In the absence of interaction F = 1 and the divergence in the integral ([?]) is cut at times t ~ th- In this case from 
((?]) we reproduce the well known result 

5a{ 0) = (8) 

For large th (i.e. for H — ► 0) the result rtq) diverges and the effect of interaction should be taken into account. 
Provided in the long time limit the function F decays faster than \j \ft the integral ^ converges even for 1 /th = 
and we get 



5ai = -a—^/DTp, (9) 
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where the prefactor a ~ 1 which depends on the function F. The precise definition of a is of little practical interest 
since this prefactor can always be removed by rescaling of t v . Of importance, however, is to describe the behavior 
of the function F(t/r v ) at t <~ t v . This allows to determine the magnitude of the dephasing time 7%. Clearly, a 
nonperturbative analysis in the interaction is needed in order to determine the function F at times t ~ t v simply 
because there exists no small parameter in the problem. E.g. if one would formally decrease the interaction strength, 
the magnitude of the dephasing time t v would increase, but one would never avoid the necessity to determine the 
function F at t ~ r v . Thus the problem of finding the decoherence time in disordered conductors is nonperturbative 
for any interaction strength. In this respect the statement of Ref. [It] that "since t v is much longer than the elastic 
scattering time, the dephasing is weak and there is no need to invoke nonperturbative ideas" remains puzzling to us. 
Indeed, at times of order t v the dephasing is strong by definition, and it is not clear how the condition r v 3> r e might 
help to avoid "nonperturbative idjsas" . 

Observing this problem AAGEj suggested to consider the limit of strong magnetic fields th <C t v , for which the 
integral (^) converges already at times t ~ r# much shorter than t v . In this case the weak localization correction 
can be calculated pcrturbatively in the interaction or, equivalently, by means a short time expansion of the function 
F. The recipe to evaluate the dephasing time from the perturbation theory suggested by AAG can be summarized 
as follows. 

In the zero order in the interaction we have F = 1 and the magnitude of the weak localization correction (^) 
increases as ^Jth with increasing th- If, expanding in the interaction, one would recover the term oc s/th, this term 



4 



could just be added to (||) and interpreted as an interaction-induced renormalization effect of the bare parameters. 
The presence of such a term would imply that F is not anymore equal to one but acquires some interaction correction. 
Nevertheless no time dependence of F and, hence, no dephasing occurs in this case and therefore the terms oc y/fn 
are not "dangerous" . If, however, the first order conductance correction is found to increase with th faster than yJf~H 

and to have an opposite with respect to Sa^ (||) (i.e. positive) sign, this would already mean that the function F 
depends on time (decays with increasing t) due to interaction and, hence, nonzero dephasing occurs. Then, if such 
"dephasing" terms are recovered within this perturbative procedure, one should look at a temperature dependence 
of such terms. If these terms are present at a finite T, but decrease and vanish as temperature approaches zero this 
would imply that interaction does not cause any dephasing at T = 0. If T- independent positive terms growing faster 
than -y/Tu are recovered one would be able to conclude that nonzero dephasing occurs at T = already within the 
first order perturbation theory in the interaction. 

We are going to demonstrate that the above perturbative strategy in principle cannot be used to correctly obtain 
the dephasing time t v for any magnetic field even though the correction San can be evaluated perturbatively in the 
limit t h <C t v . 

To begin with, we note that already the terms oc ^Jth can easily cause troubles provided they give a (positive) 
contribution to 8<j\(H) large as compared to the magnitude of the zero order term (||). In fact, the presence of terms 
oc yjTn just implies that their time dependence saturates already at short times t < th- If this saturated value 
turns out to exceed the zero order term, this would only indicate the breakdown of the perturbation expansion in the 
interaction and, hence, no definite conclusion from this expansion can be drawn. 

An even much more important problem is that the form of the function F(t) in (^) cannot be recovered from 
the perturbation theory at all. It is quite obvious that the first order perturbative terms will depend only on the 
derivative F'(0). Although in the limit th t v the value Sai(H) can be calculated perturbatively in the interaction, 
this would yield no information about the dephasing time t v . Such information can be extracted only if one assumes 
some particular form of the function Fit/r^). But this form should be found rather than assumed. This task can be 
accomplished only if one goes beyond the perturbation theory. 

Let us consider several different functions F(t/r v ). Perhaps the most frequent choice of this function is based on 
the assumption about purely exponential decay of the phase correlations, in which case one has 

F(t/T ip )=exp(-t/T v ). (10) 

As it was already discussed above, this form of the function F follows directly from the golden rule approximation. 
Substituting ( jToj ) into ^ in the limit of weak magnetic fields th 3> t v one immediately arrives at the result for the 
weak localization correction of the form (^) with tjj substituted by (1/r^f + 1/t v ) _1 . In the limit of weak magnetic 
fields th ^> t v the result (||) with a = 1 is recovered. In the opposite limit th <C t v eqs. (Q) and ( [To| ) yield 

(0) e 2 VDt 3 J 2 

i-K T v 

where Saf 1 ^ is defined in (||) . Another possible choice of the function F can be 

F(t/T v )=exp(-(t/T v ) 3 / 2 ). (12) 

The reason for such a choice will become clear later. The substitution of ( |l2| ) into (Q) again yields the result @ (with 
a = 2r(l/3)/3y7r ~ 1.0076, r(x) is the Euler gamma-function) in the limit th 2> t v , while in the opposite limit 
th *C t v from (Q) and ( |l2| ) one obtains 

3/2 



, , (o) e y/DT H ( t h 

Comparing ( |i"T| ) and (|l3|) we observe that for strong magnetic fields the interaction corrections to the leading order 
term (||) are different depending on the choice of the function F, even though for weak magnetic fields both choices 
© and (H|) yield the same result (^) with only slightly different -values of a numerical prefactor a. 
The magnetoresistance data are frequently fitted to the formulatl 

■Ax(t v /th) 



8a 1 {H) = -jDT ~: ^ »> (14) 

7T Al {T V /T H ) 

where Ai(x) is the Airy function. In the limit th S> t v this equation again reduces to (Q) with the factor a = 
— Ai(0)/Ai'(0) ~ 1.372. In the opposite limit th -C t v one finds 
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2 m — / \ 3 / 2 
(o) _ e v Dth ( t h x 



We observe the equivalence between ( |13[) and (|15| ) up to a numerical prefactor of order one. 
Finally, let choose the trial function F in the following form: 



F(i/7v,) = -=J=^=, (16) 
^(l-e- 6 */^) 



where 6 is a numerical coefficient of order one. Combining ([?]) and ([!(]) we find 



i r 



fa x (g) = ^ ^ \ \ ■ (17) 

* ^r(i + i(i + ^)) 

In the limit th <C t v one can expand this equation in powers of th / t v and get 



(0) 4-6 e 2 ^Dr H t h 
foi -M * — — (18) 

while in the limit 7~h S> from (R), ( |l6| ) one again recovers eq. (^J) with a slightly modified numerical prefactor 
a (which now also depends on the value 6). Absorbing a by a proper redefinition of t v one can plot the result (pi) 
for the trial functions (|l0|), (|l^) (for different values of 6) and eq. ( [Ti]) depending on the magnetic field (or th). 
These plots are presented in Fig. 1. We observe that all four plotted functions are very close to each other (e.g. the 
maximum deviation between jl^ ) and Sa\(H) obtained from ( |l6| ) with 6 = 4 does not exceed 0.01). If one would 
fit the experimental data for the magnetoconductance with any of these four functions one (i) would not be able to 
distinguish between them within typical error bars and (ii) would obtain the same value t v for all these functions (up 
to a prefactor a ~ 1 absorbed in t v anyway). In other words, the results for t v extracted from fitting the experimental 
data to several different functions F(t/r v ) will be practically insensitive to the particular form of F as long as its 
decay at long times is sufficiently fast to provide an effective cutoff for the integral (0) at t ~ t v . 

At the same time if one tries to extract t v from the perturbation theory in the interaction one immediately arrives 
at ambiguous and contradictory results. Let us, for example, consider the perturbative result of AAG 

(5(71 - Sa[ 0) cx Tt% (19) 



(see e.g. eq. (4.8b) of Rcf. [L3|) and, following the above paper, assume an exponential decay of correlations (|1C|). In 
this case the dephasing time t v is obtained from a direct comparison of eq. (pd| ) (equivalent to eq. (4.3b) of Ref. [l3| 
or eq. (3) of Ref. [□]) with eq. @ (or eq. (4.8b) in Ref. |l|). One obtains 

1/t£ ag « (20) 



(cf. eq. (4.9b) of Ref. 13). The result (|20|) is essentially based on the assumption about a purely exponential decay 
(|l0|). Note, however, that a-priori there is no reason to assume such a decay. [Just on the contrary, it will be 
demonstrated below that this is not the case for the problem in question.] The cutoff functions (|l^) , ( |l6| ) (and many 
others) yield the same result (||) as the function |l(]) and one can hardly make a distinction between them from the 
magnetoconductance measurements (Fig. 1). 

For instance, if one sticks to the choice (^2|), one should extract t v by comparing eqs. (|l3|) and (|l^). This comparison 
yields t v independent of th andpil/r^ cx T 2 / 3 . The latter form coincides with the well known result by Altshuler, 
Aronov and Khmelnitskii (AAK)H but is in an obvious disagreement with (|20|). If, instead of (|l2|), one uses the trial 
function (16) and compares ( |T^ ) and (p0|), one finds t v cx (4 — 6)/(Tr^/ 2 ), i.e. positive, zero and even negative (!) 
dephasing times respectively for 6 < 4, 6 = 4 and 6 > 4. However, all these dramatic differences in the first order 
results have no effect both on the form of the magnetoconductance Sad(H) (Fig. 1) and on the value t v extracted 
from it. 

Thus, whatever result is obtained in the first order perturbation theory in the interaction, it is yet insufficient to 
draw any definite conclusion about the dephasing time t v . The problem ia. essentially nonperturbative and should be 
treated as such. The corresponding analysis was developed in our papem and will be extended further in the next 
section. 
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FIG. 1. Magnetoconductance curves corresponding to different cutoff procedures: (a) purely exponential cutoff (|10|); (b) 
Airy function of eq. (0)); (c) the cutoff function F is given by the eq. ( [l6l) with b = 4 (a ~ 1.479) corresponding to "zero 
dephasing" in the first order perturbation theory; (d) the cutoff function (ftm with b = 5 (a ~ 1.582) which yields "negative 
dephasing" in the first order perturbation theory. Here t v is rescaled to absorb the factor a in eq. ([]). 



III. WEAK LOCALIZATION CORRECTION: NONPERTURBATIVE RESULTS 

In order to provide a complete description of the electron-electron interaction effect on the weak localization 
correction (||) it is in general necessary to calculate both the function fd(t) in the exponent of @ and the pre- 
exponential function Ad(t). An important observation is, however, that information about the effect of interaction 
on Ad(t) is not needed to correctly evaluate the dephasing time t v . It suffices to find only the function fd(t) which 
describes the decay of correlations in time and provides an effective cutoff for the integral (](]) at t ~ t v . The role of 
the pre-exponent is merely to establish an exact numerical prefactor. Since is defined up to a numerical prefactor 
of order one anyway, it is clear that only the function dd{t) - and not Ad(t) - is really important. 

In the subsection A we extend our previous analysisa and evaluate of the function fd(t) keeping all the subleading 
terms. This procedure is important ip-at least two aspects: (i) it allows to unambiguously settle the issue of unphysical 
divergences which was argued by VAE3 to be a problem in our previous calculatiorfl and (ii) it is necessary to establish 
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a direct relation between our nonperturbative approach and the perturbation theory in the interaction. In the 
subsection B we will perform a semi-quantitative analysis of the effect of interactions on the pre-exponential function 
Ad(t). In the subsection C we will demonstrate that the short time perturbative expansion of both the exponent and 
the pre-exponent at low T not only can (section 2B) but does lead to missing an important information about t v in 
disordered conductors. 



A. Exponent 



The function fd(t) can be evaluated by means of the path integral formalism!. This procedure amounts to calculating 
the path integral for the kernel of the evolution operator 

J ~ [vrfvp exp{iS - iS' - iS R - S7) (21) 



within the saddle point approximation on pairs of time reversed paths and to averaging over diffusive trajectories. 
Here So and S' represent the electron action on the two parts of the Keldysh contour, while iSu + Si accounts for the 
interaction. The effective action ( ^ll) was d erived in our Ref. [| for the sake of convenience we reproduce the explicit 
expressions in Appendix A (eqs. ( A16- A1S| ) ) together with the expression for the c ondu ctance of a disordered metal 



in terms of the kernel of the evolution operator J and the electron density matrix ( |A15| ). 

The saddle point approximation procedure was d escribed in details in Ref. 0. One can demonstrate that the 



contribution of the real part Sr of the action (Al£) vanishes on any pair of time reversed diffusive paths. By no 
means this cancellation occurs by chance, rather it is a generic property of a wide class of influence functionals 
describing dissipative environments. E.g. similar cancellation is observed in the Caldeira-Leggett modelo, the relation 
to which will be discussed in Section 5. Also we would like to point out that for any pair of paths Sr contributes 
only to the real part of the effective action and, hence, can never cancel an imaginary contribution from the ter m S r 
( |A19| ). Anyway, the function fd(t) in the exponent is determined solely by the imaginary part of the action Sj ( |A19| ) 
and it is given by the following expressions: 

t t 

fd(t) = e 2 / dti ( dt 2 (Ifa - t 2 , r(ti) - r(t 2 j) - I{t x +t 2 -t, r{h) - r(t 2 ))) diff , (22) 



where the function / is defined in (A14). In equilibrium it is expressed in terms of an imaginary part of the inverse 
effective dielectric susceptibility l/e(w, k) multiplied by coth(w/2T). The first term ( f22| ) describes the contribution of 
the self energy diagrams (diagrams (a) and (b) in Fig. 2), while the second term is due to the vertex diagrams ((c) 
and (d) in Fig. 2). In order to evaluate the function fd(t) ( p2|) we introduce the Fourier transform of the function / 
and then average over the diffusive trajectories with the aid of a standard replacements (exp[ik(r(ti) — r(i 2 ))])diff = 
exp(— Dk 2 \ti — t 2 \). After the integration over fc we obtain 

. . , Ae 2 D 1 - d ' 2 ( f d d x \ f duj duo' 

M*J = - to-\A I / ttzj ) / "ToZvT 



o- d {2ir) d \J l + x A JJ (2tt 

, Krl/O 9 , w cos ujt — cos to't 
V 7 wcoth — 



i /irf/9 9/ IS , LU — Oj' 1 — COSOjt 

\oj'\ d/2 ~ 2 {uj - w') coth - 



2T 



2T uj 



a 



(23) 



Here again the first and the second terms in the square brackets are respectively from the self energy (Fig. 2a, b) 
and the vertex (Fig. 2c, d) diagrams. For Id and 2d cases the integral of the first term over u>' diverges at u>' — > 0. 
However, it is easy to check that this divergence is exactly canceled by the second term, the whole integral is finite 
in any dimension and does not require artificial infrared cutoffs. Various divergences are rather inherent to the 
perturbation theory in the interaction and - at least in part - are due to insufficiency of the perturbative expansion 
in our problem, especially at low temperatures. It is also useful to note that at T — > the leading contribution to 
fd{t) in the long time limit is insensitive to a divergence contained in the first term in the square brackets (^3|) and 
can be derived only from this termQ. 

The integrals in (|2^) can be handled in a straightforward manner. Technically it is sometimes more convenient to 
perform calculations in the real time rather than in the frequency representation. Here we will present the calculation 
for the Id case. 



First we find the explicit expression for I(t, x) ( A14) 
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. 1 f dudk uj uj 



27T(Tl V dt 



: V[ttT coth(TrTt)] 



(24) 



Here "P stands for the principal value, i.e. V[ttT coth(7rTi)] is a distribution rather than an ordinary function divergent 
at t — ► 0. For a given diffusive trajectory (and on sufficiently long time scales) this function can be replaced by the 
following function of time 



{lit - 0, x(t) - x(0))) diff = / dx I(t, a02>(|t|,a;) = 



1 



1 Dt 



2iro\ \ dt 



\/— — — -V[ttT coth{irTt)] 

1X(T\ V 7T V dt 



-—V[TrTcoth{TrTt)} / dx 



\x\ exp 



'AD\t\ 



V^D\t\ 



Substitution of this equation in eq.(22) yields 



2 ' 



Dh 

7T 



D\2t 1 -t 2 -t\ 



TV 





Let us first integrate this expression over t\ and then integrate the result by parts. We obtain 

t 



(25) 



(26) 



hit) 



2e 2 D 



7TCT1 



— / dt' nT coth(irTt') ( — 



-Vt 7 + Vt-f 



(27) 



2T e /7T 



The short time cutoff in ( p7| ) is equivalent to a sharp cutoff at uj — 1 jr e in the frequency domain. 
In the quantum regime irTt <C I we find 



e 1 2D 2e 2 Dt f 2irt 
hit) = J 1+ \ — In 6 

TT&l V T e 7T(Tl V 7T V T e 



TxTt < 1. 



(28) 



Note, that apart from the leading linear in time term there exists a a smaller term cx y/tln(t/T e ), which also grows in 
time. 

In order to find the function f\{t) in the opposite thermal limit irTt ^S> 1, let us rewrite the integral j27j ) in the 
form 



Mi) = J— \ .r^-t - 7rTt 3/2 coth(7rTi - 1) + fV^T f — fcothx - I 

TTCTl V 7T V 2T e J y/X \ 

\ o 

TTTt t 

J dxy/x(cothx - 1) + y dt' ttT coth(nTt')Vl~~F 



sinh 2 a; 



2VttT 

Making use of the following integrals 
coth x — 1 



sinh 



2r e /7r 



2^(2/ 



efe-v/a^cotha; — 1) 



2 V 2H2 



(29) 



we get 



f /" -'--^3/2 , ^(Vg) ^ 3C(3/2) I 



7r<Ti V 7r ^ y 2r, 
/ I 



Viln 



\4Tr e 



/ 1 



V2 



4V2 







\TVt 



2irTt 3/2 e- 27rTt + C(Vte" 27rTt )| , nTt > I. 



(30) 



We observe that in both cases (|28|) and fl30| ) there exists a linear m lime temperature independent contribution to 
fi(t) which determines the dephasing time t v at low temperaturesaQ. Beside that at Tt » I there exists another 
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term cx Ti 3 / 2 which yields dominating contribution to t v at high temperatures T >T q ^ l^T v T e , where the result of 
AAkI t v oc T~ 2/3 is recovered. 

In addition to both these important contributions all four diagrams of Fig. 2 yield subleading terms in the expression 
for fi(t) which also grow with time, albeit slower than the main terms. These subleading terms also contribute to 
dephasing even at T = (cf. eq. (|2|)), however this contribution is always smaller than that of the leading terms, 
typically in the parameter yjT e /T v . This result is in contrast with the statement of Ref. [16], where it was argued that at 
t ~ t v the contribution of the vertex diagrams to f\ (t) can be comparable to that of the self-energy diagrams and the 
term oc t/^/f^, which is the most important at T — » 0, can be canceled. A straightforward calculation demonstrates 
that this is not the case. 

A similar calculation can be performed in 2d and 3d dimensions. In any dimension the result can be expressed in 
the form 

f d (t)=t/T° + 5f d (T,t), (31) 

where we defined 

1 K d e 2 (2D) 1 - d / 2 



t-0 d/2 
> TTCTdTe 



(32) 



A numerical prefactor n d in ( p2[ ) (determined for a sharp high frequency cutoff at u) = l/r e ) is K\ = 1 for Id, k<i = 1/4 
for 2d and K3 = 1/3tt for 3d. A detailed expression for the function 5fi(T,t) is given in (|2|) and (|30|), where we 
retained also several subleading terms needed for further comparison with perturbative results (Section 4). For 2d 
and 3d systems we will present only the leading order contributions to 8f d (T,t). In 2d we find 

<5/ 2 « = — Tt«l, 

<5 h{t) = ^Ttln(Ti), Tt » 1, (33) 

Here 70 ~ 0.577.. is the Euler's constant and t 3> r e . Similarly, in the 3d case we obtain 

5/ 3 (t)~7.8^ Ti<l, 
(Tt ) 3 / 2 

5f 3 (t)^3 [ e > t Tt»l, (34) 
In 3d we used a standard approximation and replaced lu coth(w/2T) — \lu\ by 2T6(T — |cij|). 



B. Pre-exponent 

As it was already discussed above, the pre-exponential function A d (t) does not play any significant role in our 
problem. Therefore its rigorous calculation at all times (which is a separate and quite complicated problem) will not 
be discussed here. Of importance is to qualitatively understand how the function A d (t) is modified in the presence of 
the electron-electron interaction. Therefore we will restrict ourselves to semi-quantitative arguments which, however, 
turn out to give surprisingly good agreement with the rigorous results obtained in Sec. 4 in the short time limit. 

It is well knownB'Ea that without interaction the function A d (t) is related to the return probability of diffusive 
trajectories to the same point after the time t. In the presence of dissipation (described by the term Sr in the 
effective action) the particle energy decreases and its diffusion slows down. This implies that at any given time t the 

function A d (t) should exceed the pre-exponent A d °\t) evaluated without interaction. On the other hand, at least if 
the interaction is sufficiently weak, diffusion will still take place at all times and, hence, A d (t) will decay in time, 

albeit somewhat slower than A d °' (t) . 

Now let us try to find a typical time scale at which the deviation of A d (t) from A d °\t) becomes of the order of 

A d °\t). For the sake of definiteness we restrict our analysis to the Id case. As we have already discussed, the real 
part of the action Sr vanishes on the time reversed diffusive paths. In order to evaluate the contribution of Sr in the 



path integral (|2lj) or (A16) we need to include fluctuations around the time reversed paths. We assume that thes e 
fluctuations are small and neglect them in the arguments of the functions R(t\ — t 2 ,ri(tj) — rk{t n )) in eq. ( A1S| ). 
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These fluctuations are, however, important and should be kept in the arguments of the functions 1 — 2n(pi,ri). In 

equilibrium one has 1 — 2n(pj , r j ) = tanh ^ , where we defined £j = + U(rj ) — /i. Within the above approximation 
we get 



S R ~— dh dt 2 {[(R(h - t 2 ,r{t 1 ) - r{t 2 ))) diS - {R(t! - t 2 , r(h) - r(f - t 2 ))) 



/diffJ 





x [tanh(&(i 2 )/2T) - tanh(6(t 2 )/2T)]} 



(35) 



In addition to the contribution (|35|) one should also take care of the corrections to the action So (A17) due to the 
interaction. These corrections turn out to be of the samg order as (|35|). In the presence of interaction the classical 
paths change and satisfy the following Langevin equatior 



mr(t') + VU(r(t')) + e 2 dt"V ' r (t')R{t' ~ t" , r(t') - r{t")) tanh 



an 

2T 



-eE(t',r(t% 



(36) 



where E(t',r) is the fluctuating electric field due to the Nyquist noise. From this equation we find 

t ti , t 

Solp,r]~Solp {a \r^]+ J dh J dt 2 le 2 J dt"r(t 2 )V r(t2) R(t 2 - t", r<°> (t 2 ) - r<°) (*")) tanh ^ 

* 

^So[p ( V (0) ]+ / dtJe 2 J dt" R(h t" , (h) - r<°> (t")) tanh ^ \ 

\ ' diff 

t t 

c S [p ( V (0) ] +e 2 / dtx / dt 2 (jZ( tl - i a)r (0) (il )_ r (0) (t2)) ^ tanh^. 



diff 



(37) 



Here r^ ^ is a classical diffusive path without interaction. The energy £ is conserved along such a path. The average 
energy change due the noise field E(t',r(t')) (|3^) vanishes and therefore was omitted in (p7|). 
Adding eqs. ( j35[) and (|37| ) together we obtain 



(38) 



(39) 



2 

So - S* - S fl ~ S \p uri ] - S [p 2 ,r 2 ] + 1- JdhJ dt 2 K [tanh(£i(i 2 )/2T) - tanh(£ 2 (i 2 )/2T)] , 



o o 



where we defined 



ft = (i?(ii - t 2 , r(ti) - r(t 2 ))) diff + (E(ti - ta, r(t x ) - r(i - i 2 ))) 



diff ■ 



Within our simple approximation the paths r± and r 2 in the action Sr are considered to be independent from each 
other. Therefore the kernel (Oil) can be expressed in the following form: 



J = U{t,r lf ,r li )U + {t,r 2f ,r 2i )e-^ t \ 



where 

U(t, ri/, r u ) = J T>r J Vpexp 
It is convenient to define the following function: 



iSa[p,r] + i— / dt± / dt 2 Ht&nh 



6(t 2 ) 

2T 



o o 



e 2 /2iJ e 2 /L>i 



(40) 



(41) 



(t) = — dh / cft 2 ft~ W — 

7r<Jl V Te 2(7 1 V 7T 



(42) 



Then the operator (^) can be rewritten as follows: 
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C/(t,r 1/ ,r li ) = X! e " i£At+ ™ (t)tanh| * MrifWxiru). (43) 

A 

Here tp\ and £>, are respectively the eigenf unctions and the energy eigenvalues of a single electron Hamiltonian 
#o = + (r) - 

In the absence of the interaction the pre-exponent is given by the following expression: 
According to eq. fli^ ) the products £ji should be replaced by £jt — u(t) tanh J§^. Thus we get 



Ai(t) = -4= / «i / «2 f J- tanh 4 ) LV " ~ (44) 



(£i — £2)* — u(t) ^tanh |^ — tanh J|? 



Let us emphasize that the estimate (Q) was obtained with the aid of several crude approximations and, in particular 
in the long time limit, corrections to this simple result can easily be expected. However, since we are not interested 
in the details of the long time behavior of Adit), the result ([5]) is already sufficient for our purposes. The main 
properties of Ad(t) are as follows. 

Firstly, eq. (0) determines a typical time at which Ai(t) significantly deviates from A[ (t). This scale (which we 
will denote as ta) is set by the function u(t) and at low T can be determined from the condition u(ta) ~ 1. Combining 
this condition with eq. (|J) and observing that the first term in this equation equals to t/r° and the second term is 
small for all times t ^> r e , we conclude that - at least for sufficiently low temperatures - the time scale ta is of the 
same order as the dephasing time at T = ((H), i.e. ta ~ t®. Thus for all t < t® the effect of the interaction on 

the pre-exponent is small and for such tunes one can safely approximate Ai(t) ~ A\ (t). This approximation was 
already used within our previous analysisQ. 

Secondly, the estimate (|44| ) illustrates again an intuitively obvious property of the pre-exponent: in the long time 
limit A\(t) decays in time. Thus no compensation of the exponential decay of correlations oc exp(— fi(t)) can be 
expected from the pre-exponcntial function A\{t) at long times. Hence, in our problem the effect of interaction on 
the pre-exponent A\(i) can be disregarded also in the long time limit t > i~. 

The same analysis can be repeated for 2d and 3d cases and the same conclusions will follow. 



C. Discussion 



Our consideration allows to suggest the following transparent picture. The dephasing time is fully determined by 
the imaginary part of the effective action Si which contains "coth" . In other words, the function in the exponent of 

©is 

fd(t) = f d [Si}. (45) 

The real part of the effective action Sr (which depends on "tanh" and contains information about the exclusion 
principle) contributes to the (unimportant for dephasing) pre-exponent Ad(t) in (^), i.e. 

A d (t) = A d [S R }. (46) 

The splitting between the exponent and the pre-exponent of the type (^), ([l6|) holds also for the exactly solvable 
Caldeira-Leggett model. This will be demonstrated in Sec. 5B. 

Although the difference between Ad(t) and A d (t) cannot have any significant impact on the dephasing time r v , 
this difference should be taken into account if one evaluates the weak localization correction perturbatively in the 
interaction. In the limit th <C t v a short time expansion of both the exponent and the pre-exponent is sufficient and 
the weak localization correction can be represented as a sum of three terms 

Sa d = 6a d 0) + S<r dcph + <5a tanh , (47) 

where Sa d °^ is the "noninteracting" correction and 
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S(J deph 



2e 2 (DT„) 1 - d / 2 

7r(47r) d / 2 



dx- 



,d/2 



fd(XT H ), 



(48) 



2 1 d/2 °° 



Making use of the result (^) for the Id case in the limit T — > and at short times we obtain 

= C ° SK ^ nKt)] - + «(t) + 0(«»(t))]. (50) 

We will keep track only of the leading contribution to the function ( ff2| ) u(t) ~ V T ^- This is sufficient within the 
accuracy of our estimate . Combining ( [48| ) and ( fi9| ) we observe that the sum of the last two terms in (^) depends 
on the combination 



A (i) -«(*)• 

The term t/r^, drops out of this combination, it is contained both in fi(t) and u(t) and cancels out exactly. The same 
cancellation occurs in 2d and 3d cases. This cancellation illustrates again the conclusion of Sec. 2: it is impossible to 
obtain correct information about the dephasing time even from the correct first order perturbative analysis. 

The accuracy of our estimate of the pre-exponent at short times t < t® can also be checked by means of a direct 
perturbative calculation. This calculation is presented in the next section. It demonstrates that the above cancellation 
of the first order linear in time T-independent terms from the exponent and the pre-exponent has a general origin 
and is not related to the quasiclassical approximation and/or disorder averaging at all. 

The final results for the weak localization correction to the conductance presented in the next section are mainly 
focused on the Id case. Here we provide the results for the 2d case. In the "perturbative" limit th <C t v one obtains 
from @, ©-© 

6a 2 (H) - 5af\H) = ^^Tr H \n(Tr H ) (51) 
in the thermal limit Tth 1 and 

in the quantum limit Tth-j=C 1. Here Rn is the sheet resistance of a two-dimensional film. The result ( |5l| ) coincides 
with that found by AAGta in the limit Tth ^> 1. An opposite limit of low temperatures was not considered in Ref. 
O at all. We will perform a detailed comparison of our results with those of AAG in the next section. 



IV. PERTURBATION EXPANSION 



Now let us analyze the expression for the weak localization correction to the conductivity bad perturbatively in the 
interaction. The structure of this section is as follows. In the subsection A we will derive general exact results for 
the system conductance in the first order in the interaction. In the subsection B we will demonstrate that the exact 
first order diagrams do not cancel at T — and, moreover, that the result cannot in general be interpreted as an 
effective renormalization. We will also demonstrate that some previous statements about an exact cancellation of the 
first order diagrams in the limit T — > are nothing but artifacts of insufficient approximations, the main of which 
is the golden rule approximation. A detailed calculation of the weak localization correction in the first order in the 
interaction is performed in the subsection C. There we will also identify the contributions to this correction coming 
from the exponent and the pre-exponent, see section 3. In the subsection D we will present a detailed comparison of 
our analysis with that developed by AAG in Ref. [l3|. 
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(a) 




(c) 




Jp j« t. 



(b) 




(d) 



FIG. 2. The four first order diagrams. The time ordering is as follows: ti > ti > ts > 0. 



A. General Results 



The p erturbation theory can be constructed by means of a regular expansion of the kernel of the evolution operator 
J ( Al£) in powers of iSr + Si. In the first order one obtains four diagrams presented in Fig. 2. The contribution of the 
self-energy diagrams (Fig. 2a, b) was analyzed in details in Appendix A. It was demonstrated that this contribution to 
8a can be written in terms of the the evolution operator U for noninteracting electrons. The corresponding expression 
is defined in ( A12- |A14 ). It is equivalent to the result ( |Al[ ) obtained diagrammatically by AAG. 

Let us express the evolution operator U in the basis of the exact wave functions of noninteracting electrons: 



U(t,n,r 2 ) = (ri\exp(-iHot)\r 2 } = J2e- i ^ t tPx(riW x {r 2 ) 



(53) 



Obviously the representation (|53j) holds both with and without the external magnetic field with the only difference 
that in the latter case the energy levels are doubly degenerate, while in the former case this degeneracy is lifted by 
the magnetic field. 

The density matrix p which enters the expression (A12) can also be expanded in the basis of the eigenfunctions ip^. 
We find 



(1 - 2p)(n,r 2 ) - ^tanh^^A(ri)^(r 2 ), 

A 



(54) 
(55) 



Now let us substitute (|53|-p5|) into the expression ( A12 ). Performing the two time integrals after a straightforward 
algebra (see Appendix B for details) we obtain the correction to the conductivity due to the self-energy diagrams of 
Fig. 2a,b 



5a. 



a/3 



— / *i y 

2Va J /-J 

n A1..A4 



2tt 





coth 



tanh — - 



-^+tanh^ 
2T 2T 



( 1A4A1 -AiA 2 , -A4Ai 'AiA2 1 



M 



A2 A3 ; A3 A4 



(56) 
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where we defined the matrix elements 

i^ A = = J drr Xl ( r )j^x 2 (r), M x * x *** x * = 
and the function 



rl 3 k 1 

' ' T7 (A 2 |e^|A 3 )(A 3 | e -^|A4) 



(2tt) 3 k 2 



F(ti,w,f A i»6u)= / *2 /" <ftscos(-fAi*i+C)ta*3 + (U+ w )( t 2-*a)+64(*i-*a)) 



(57) 



(58) 



o o 



COS £21 i 



£42[cOs((£ 32 + U)ti) - 1] + (f 32 +w)[l - COS £ 42 tl 



sin £ 2 i*i 



£42 (£42 - 62 - ^)(^32 + w) 

£42 sin((£ 32 + u)ti) - (£ 32 + w) sin£ 42 £i 



^42 (^42 - 62 ~ ^)(C32 + W) 



Here we introduced the notation £jj = £a 4 — £\, . 

The term Sa^g in ( [56| ) describes the correction due to the non-screened Coulomb interaction. It is defined by the 
following expression: 



+00 



-|v J 1 dt * E 

n Ai,A2,A4 



+ ^9 Jc 



x (A 2 | [1 , gfelg) |A 4 ) sinh^ti + £a 2 * 2 + £a 4 (*i ~ *a)]. 



(59) 



\n — T2\ 

The contribution to 5a from the vertex diagrams of Fig. 2c, d can be found analogously (see Appendix B). We get 



X vert 



OO 

2Va J X X 

n A1..A4 





du> 
2^ 



tanh — - 

d£ Al 2T 



X2X3 -XiXi 1 -A 2 A3 -AiA4\ 
3/3 +3 13 3 a j 



coth — + tanh — 
2T 2T 



M 



A3 Ai ; A4 A2 



(60) 



Here we have introduced the following function: 



ti 



G(*i,W,& 1 ..&J = J dt 2 J dt 3 COS(-£x 3 (tl-t3)-£\it3+£\j2+£\ 2 (tl-t2)+w(t 2 -t 3 )) 


(?42 + £3l) COs[(£ 3i - £ 42 ~ W)*l] - (£l2 + w) C0s£ 42 ii - (£31 ~ w) COS £31^ 



(61) 



(62 + £31) (£31 - w)(&L2 + w) 
(£42 + £31) sin [(£31 - £42 - u)h] + (£42 + u) sin£ 42 ti - (£ 3i - uj) sin£ 3 i ti 
(£42 + £31 ) (£31 - w) (£42 + oj) 



cos £ 2 i<i- 



-sin£ 2 ii 



Despite an obvious similarity in the structure of the self-energy (eq. (p6[)) and the vertex (eq. (|60|)) corrections to the 
conductivity these two expressions differ in several aspects: the terms containing coth(w/2T) in (^fjj) and (|60| ) have 
the opposite signs, the functions F ( |58| ) and G ( |6l| ) of the energy arguments £a are different and the matrix elements 
entering (|5^) and (^) depend on different indices. 

It is important to emphasize that the eqs. (|56]-|6l|) determine the total correction to the conductivi ty te nsor 
&af3 = + ^ap" wn i cn is identical to the initial results ( |Al[ ) and ( A12 ). In deriving (|56]-|6l|) from ( A12 ) no 
quasiclassical approximation, no averaging over disorder and/or no other approximation of any kind has been made: 
the above equations are exact quantum mechanical results in the first order in the interaction. Therefore these 
equations can be conveniently used to test the statement about_the full cancellation of .the .first order diagrams at 
T — ► which is quite frequently made in the literature (see e.gBj as well as recent worksEJiia and further references 
therein). 



15 



B. Breakdown of the Fermi Golden Rule Approximation 



Although here we are mainly interested in the contribution of the diagrams of Fig. 2 to the current-current 
correlation function, the structure of the result is by no means specific to this function only. The very same structure 
- perhaps aparjt-JMpi, the matrix elements of the current operator - is reproduced if one calculates e.g. the inelastic 
scattering timeE3'E3ll3 and similar quantities. This is q uite natural because the results for different quantities follow 
from the expansion of the same evolution operator J (AL6) in the interaction. Hence, the analysis to be presented 
below is general and can be applied to various physical quantities evaluated by means of the diagrams of Fig. 2. 

Let us consider the self-energy diagrams of Fig. 2a f bj Just for the sake of clarity let us repeat the statement we are 
going to test: according to Fukuyama and Abrahamsc3 and to some other authors the contribution of these diagrams 
vanishes in the limit T — > because the result contains the combination 



d , e \ 

& tanh 2r) 



, uj e — uj 

coth — - + tanh — — — 
2T 2T 



(62) 



under the integrals over e and ui. This combination restricts both integrals to the regions |e| < T and 
makes the result to vanish completely at T — 0. 

Already the first inspection of the expression (|56],p8|) allows to observe that it is the combination 



< T and 



^-tanh^ 



coth h tanh 

2T 2T 



(63) 



and not (62) which enters the exact quantum mechanical result. This combination is not zero even at T — because 
£\ 3 7^ £ai — high frequencies \u>\ > T do contribute to the integral and, moreover, this integral may - depending on 
the spectrum of the fluctuation propagator - even diverge for large u> unless one introduces an effective high frequency 
cutoff. We would like to emphasize that these conclusions are general and do not depend on any particular form 
of the matrix elements ( |57| ) . Thus the statement of the above papers that the contribution of the diagrams of Fig. 
2a,b vanishes in equilibrium at T = is proven to be incorrect. Below we will demonstrate that this poorly justified 
statement is a result of several rough approximations, the main of which is the golden rule approximation. This 
approximation may sometimes yield correct leading order results in the high temperature limit, but it breaks down 
at sufficiently low T. 

In order to illustrate this point let us first make a simplifying assumption. Namely, let us for a moment restrict 
our attention only to the contribution of the terms with £\ x = £a 2 = £a 4 - Below we will see that this assumption is 
not sufficient to properly evaluate the first order perturbation correction to the conductivity: in order to do that it 
is important to allow for a (possibly small) difference between £\i and £a 2 • But such an approximation is sufficient 
for calculation of some other physical quantities, like the inelastic scattering time, and we will adopt it for a moment 
just in order to demonstrate the failure of the golden-rule-type perturbation theory in the interaction. 

The contribution of the terms with £ai = £a 2 — £a 4 to the conductivity Sa sc reads 



e 2 T ?Ai= ^ =Ca4 / d £ A 
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Let us first evaluate this expression within the Fermi golden rule approximation: 

1 - cos((£ 3i +u)ti) 



fei+w) 2 



(64) 



(65) 



Substituting ( |65| ) into (|64| ) we obtain 
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(66) 
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This expression implies that within the Fermi golden rule approximation the self-energy diagrams of Fig. 2a,b yield a 
linear in time decay of the initial quantum state with the corresponding relaxation rate proportional to T . Obviously, 
such relaxation rate vanishes at T = in agreement with Ref. ^ and others. 

Now let us carry out an exact frequency integration in (E4J) without making the golden rule approximation (|65|). ft 
is fairly obvious that the integral over u) is not restricted to \u>\ < T and even diverges at high frequencies. As before, 
in order to cure this divergence we introduce the high frequency cutoff uj c l/r e . For simplicity we also assume that 
the energy difference £31 is smaller than l/r e . Then in the limit T — > we find 
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coth • 
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(1 — coswti) . 



(67) 



The first and the third terms in the second line of this expression come from coth(w/2T) while the second term 
originates from tanh(^ 3 /2T). We observe that the first two terms are the same as in the golden rule approximation 
(pq). These terms enter with the opposite signs and exactly cancel each other at T = because in this limit 
(d/d£\ 1 )taiih(£\ 1 /2T) reduces to a <5-function and therefore £a 3 = £31. The last term does not vanish even at zero 
temperature, this term is not small and obviously contains the contribution of all frequencies up to l/r e . The integral 
over u> contained in this term can be easily evaluated. We will do it a bit later when we fix the dependence of the 
matrix elements M A2 ' A3;A3 ' A4 on energies. Now it is only important for us to demonstrate that the last term in ( |67| ) 
is completely missing within the golden rule approach employed in Refs. |2(],[l7],[l8| and others. It is obvious, therefore, 
that this approachri-ails to correctly describe the system behavior at sufficiently low temperatures. 

Note, that AAGJ13 also did not observe an exact cancellation of diagrams of the first order perturbation theory in 
the interaction. However, they argued that the remaining terms provide the so-called interaction correction to the 
conductance which can be viewed as an effective (temperature dependent) renormalization of the bare parameters 
and has nothing to do with dephasing. Already from the form of the third term in the right hand side of (|67]) one can 
conclude that in general this is not true. Indeed, if one adopts that for £ai — £a 2 = £a 4 the dependence of the matrix 
elements M A2,A3;A3,A4 on the energy difference £31 has the form 



M A 2 ,A 3 ;A 3 ,A 4 ^ |£ 31 |<*/2- 



(68) 



(cf. eq. (2.33) of Ref. |l||), and integrates the product of M A2 ' A3;A3,A4 and the last term in ( |67| ) over the energy £31 one 
immediately observes that after the cancellation of the unphysical divergence (which is also contained in the vertex 
diagrams of Fig. 2c, d and enters with the opposite sign, see also Section 3a) one obtains the contribution oc ^/th\(t / T e ) 
in Id and ln(t/r e ) in 2d. This contribution is just a part of the function 6fd(T,t) (31) at T — > 0. It grows with time, 
contributes to dephasing and obviously cannot be reduced to the renormalization of the initial parameters which 
would be provided by a time-independent term. 

In order to understand why AAG arrived at such a conclusion it is appropriate to highlight the approximation 
employed in Ref. [l^. As a first step they split the total contribution to 8a into two parts, eqs. (5.12c) and (5.12d), 
effectively rewriting the combination (p3|) in the following equivalent form: 



J- tanh % 
d^ 2T 



coth — - + tanh — 
2T 2T 



tanh 



£a 3 
2T 



tanh 



2T 



(69) 



The first Affo terms in the square brackets of (69) were interpreted by AAG as a "dephasing'.!— (Contribution (eq. 
(5.12c) oJl3) while the last two terms are meant to be the "interaction" correction (eq. (5.12d) ofEJ). Obviously, the 
contribution of the first two terms vanishes at T — > 0. In order to understand the behavior of the remaining terms we 
make use of (p4), (68) and observe that the contribution of the last two terms in (|6^) is proportional to the following 
integral 



du d£\ 1 / d£x 3 



tanh — tanh ■ ^ 



2T 



2T 



1 - CQS((£ 3 1 + m)t) 

(61 +^) 2 ' 



(70) 



The approximation employed by AAG while evaluating such a combination is equivalent to ignoring the oscillating 
cos-term in (70). After dropping this term and making the integral dimensionless one can easily observe that the 



remaining integral has the form 
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Ai/VT (71) 

in Id and A2I11T in 2d, where A\£ are temperature- and time-independent constants. AAG interpreted these 
contributions as an effective renormalization due to interaction. Note, however, that it is correct to drop the cos-term 
only at sufficiently long times Tt ^> 1, while at smaller Tt < 1 this term is important. Evaluation of the integral j70| ) 
in the latter limit yields 

BtVilnit/Te) (72) 

in Id and B% ln(t/r e ) in 2d, where B12 are again temperature- and time-independent constants. It is fairly obvious 
that the term d73) already cannot be interpreted as a renormalization effect from an effective static potential. This 
term explicitly depends on time and actually contributes to dephasing! 

Now we are aware of the behavior of the integra l (f70[) at all times: at t = this integral is obviously zero, it grows 
with time as ([72]) for Tt < 1, reaches the value (|71[) and saturates in the long time limit Tt ^> 1. Clearly, in the 
interesting limit T — ► the behavior (|7l|) can never be realized, the term (f7(j|) grows at all times and contributes to 
dephasing. In this limit we are back to the result (^37|). The perturbation theory strongly diverges in this case. It 
also diverges at finite temperatures, thus the corresponding expressions can only make sense if one introduces a cutoff 
at times much smaller than the dephasing time t v , because at times t ~ t v all orders of the perturbation theory 
should be taken into account. In Ref. [l3] this cutoff time was chosen to be the magnetic-field-induced decoherence 



time th "C t v . Thus the approximation leading to the time- independent term (71) is valid only for Tth ^> 1, in 
which case the contribution ([72]) is anyway much smaller than that from the first two terms in (|69|) and, hence, can be 
safely ignored in the above limit. On the other hand, in the most interesting limit Ttr ^ 1 (which is compatible with 
Tt v ^> 1) the contribution (f70| ) dominates, its behavior is given by eq. ((72]) rather than by eq. (|7l]) and, consequently, 
nonzero low temperature dephasing is observed already in the first order perturbation theory in the interaction. We 
will come back to this discussion in Sec. 4.D and in Appendix C. 

Let us emphasize again that no approximation was done during our derivation presented in Sec. 4A. Our main 
goal here was to demonstrate that the absence of the cancellation of diagrams in the first order perturbation theory 
has nothing to do with the quasiclassical approximation and/or disorder average as it is sometimes speculated in the 
literature. 



C. Perturbative weak localization correction 



Now let us perform a systematic evaluation of the exact expressions (p6|)-(|6l|) obtained within the first order 
perturbation theory in the interaction. Our calculation consists of several steps. First we notice that the expressions 
(]56|)-(pT|) contain the full information about contributions from all energy states. Since here we are interested only in 
the weak localization correction to the conductance we should restrict our attention to the time reversed energy states 
and evaluate the matrix elements for such states. The matrix elements for the current operator can be extracted from 
the expression for the weak localization correction without interaction §cr™a. Starting from the standard expression 
for this correction (see e.g. Ref. 13) and rewriting it in terms of the matrix elements for the current (^) we obtain 
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(73) 



The expression for the matrix elements of the currents jj (£21) is readily established by comparison ([73|) to the well 
known quasiclassical result in the absence of interaction (pi). We find 



A1A2JA2A1 



4e 2 DV 



DQ 2 + 1/th 



2e A V 



.£v(l n m_ii n( i + (6lTff) 2 ; 



d d Q _ 
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Id 

2d 



(74) 



As it was already pointed out above, these expressions are only valid for the time reversed states Ai and A 2 relevant 
for the weak localization correction. For later purposes let us also rewrite the above result in the Id case in the real 
time representation: 
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rf6iii(6i)cos6i*i = — ^ — \ ?p=^- ( 75 ) 



The next step in our calculation is to identify the contribution to 8a responsible for dephasing. As it was already 
demonstrated above within the nonperturbative analysis, this contribution is determined by the function fd{t) (p3| ) in 
the exponent (^) . Clearly, in the first order in the interaction this contribution is obtained by expanding the exponent 
in @ up to the linear term in fd(t) and ignoring the effect of interaction on the pre-exponential function Adit). 
Hence, this "dephasing" contribution should have the form 

+00 

S <T = i / *i /-(*o E (^ tanh If) (^ lA2 ^t Al +^ AlA2 ^ A2Al ) ( 76 ) 

We observe that this expression contains the function cos £21^1 and does not contain sin £21^1 ■ Furthermore, from 
the above analysis we know that the function fd(t) contains only coth(o>/2T) and does not depend on tanh(^A 3 /2T). 



Therefore in the general result for the conductance correction (56)-(pl[) we will hrst take care of all terms which 
contain the product coth(w/2T) cos^i^i leaving all the remaining terms for further consideration. 

Consider the "cothxcos" terms originating from the self-energy diagrams of Fig. 2a, b. For such diagrams one 
should put £42 = 0, then from (|56|), ( j58|) one will immediately observe that the contribution of the "cothxcos" terms 
can indeed be represented in the form (uw where 



f d (t,\ 2 )-—l^J—M o;coth— fe2+w) 2 • (77) 

A3 

Let us replace the summation over A3 by the integration over £x 3 . Assuming that the matrix elements M A2A3;A3A2 
depend only on the energy difference £32, making a shift £32 + u> — > u> and denoting £32 — ► to' we find 

,.«,„/ x 47re 2 f dwduj' r , ... , uj — ui' 1 — coswi 

fT(t) = / —r w-w' coth— - . (78) 

This expression does not depend on A2 and exactly coincides with the first term in eq. ( p3|) if we identify the matrix 
element M(w') as: 

Note, that the energy dependence of the matrix elements M oc \uj'\ d l 2 ~ 2 ( |79| ) determined within the above procedure 
is in the agreement with the conjecture (|6|) as well as with eq. (2.33) of Ref. |l3| . 

The contribution of the "cothxcos" terms contained in the vertex diagrams of Fig. 2c, d can be evaluated analo- 
gously. Again one should consider only the part of the function G ( |6l| ) which contains cos(^2i^i)- For the contribution 
of the time reversed states to the vertex diagrams one should identify £42 = £31. Making use of this equation, from 



the corresponding terms in ( |60| ) , ( |6l| ) one finds 

+00 

^dcph,v Cr t = -±- j dhfr^Xr) £ (-|-tanh|0 (j^j^ + 3^3^) cos^i, (80) 

A1A2 



where 



/r rt (*.Ax) = -^£ J |^ A ^ A Woth 



UJ COSLUt — COS^3ii 

2T eiw^ ■ 



(81) 



By comparing eq. (^TJ) with the second term of the expression ( p3| ) we observe that they coincide provided one denotes 
£31 = £42 = and again assumes that the matrix elements depend only on the energy difference M A3Ai;A<lA2 = M(u/), 
where M(u/) is defined in eq. (|79|). Furthermore, in order to identify eqs. ( |76| ) (with fd(t) — > /J ert (*)) and ( j8C| ) we 
have to assume that j X2X3 j X p lXl + ^ 2A3 j AlA4 = ja lX2 jp 2Xl +3p lX2 Ja 2Xl - This completes the analysis of the "cothxcos" - 
contribution from the vertex diagrams of Fig. 2c, d. 
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Thus, we have explicitly demonstrated that the perturbative "dephasing" contribution to the conductance obtained 
before from the nonperturbative analysis can also be identified in the first order perturbative expansion provided one 
infers the matrix elements M in the form ([79]) . With this in mind one can immediately write down the expression for 
<5<7^T h in the form (fl8j). For a quasi-ld case at low T we find 
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In the opposite high temperature limit we obtain 
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(82) 



7tTt# > 1 



(83) 



Now let us come to the final step of our calculation and evaluate the remaining terms in the general result (pq)- 
(|6l|). We notice that the contribution of all terms containing the combination coth(o;/2T) sin £21*1 vanish after the 
integration over the energy £21- The same is true for the terms containing tanh(£x 3 /2T) in the contribution of the 
vertex diagrams (60). These observations imply that all the remaining nonvanishing terms come from the self-energy 
diagrams of Fig. 2a, b and contain 1 — 2p or tanh. We will denote their total contribution as £)cr tanh . We already know 
from the above analysis that this contribution comes from the expansion of the pre-exponent A\ (t) to the first order 
in the interaction. Collecting all such terms from (]56|), fl5Sj ) and J59|), we obtain 
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(84) 



where 
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As before, the above equations were obtained from the exact ones by imposing £ 42 = 0. In order to establish a 
somewhat closer relation to the approach developed by AAG we also note, that it is the contribution Sa tanh 
which contains the so-called Hikami boxes within the diagrammatic analysis of Ref. [l^. AAG argued that (partial) 
cancellation of the first order diagrams can only be observed if one takes the Hikami boxes into account. Below we 
will demonstrate that this is not the case. Actually we have already shown in Sec. 4B (eqs. (^), (|67|)) that this 
cancellation (of the linear in time "golden rule" terms only!) in the first order at T — occurs already before disorder 
averaging and thus has nothing to do with the Hikami boxes. Now we will illustrate this fact again by means of a 
direct calculation. 

Let us first consider the term Sal^} h (|8q). The integral over uj can be evaluated exactly and we get 
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Further calculation will be performed for a quasi-ld case. We also make use of the real time representation of our 
integrals, as it was already done before. For Id systems we obtain from ( |79|) 
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Also we will use the following relation 
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Substituting fl75J), (|S9|) and (|9C|) into eq.(g§), we find 
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After simple algebra this equation can be converted into the following integral: 
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In the quantum limit ttTth < 1 we get 
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To consider the opposite thermal limit itTth 3> 1 it is convenient to rewrite this equation in the following form: 
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This equation yields 
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Here we have used the following integrals: 
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Now we turn to the correction (|S6|). To begin with, we should handle a divergence which appears in the integral 
over u) for the term linear in t±. It is easy to demonstrate, however, that this divergence is fictitious. It disappears 
completely if a more accurate expression for the matrix elements is used. This expression reads: 
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Now we can use the analytical properties of the function l/e and write 
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Substitution this identity into eq. ( |86| ) we immediately observe that the term containing —1 is exactly canceled by 
the correction 5a c . Thus the result is finite and has the form 



6a c + 5al^ h = Sa\ anh + &j* anh , 
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The contribution fl98| ) can be transformed and evaluated analogously to the term Sal^ h . We find 
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After simple transformations we obtain 
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In order to evaluate the term Serif we use the following expression 
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It can be obtained by comparing the two expressions for the matrix element M(£3 2 ). The first expression, 



follows directly from the eqs. ( |22|j25| ) and (77), and the second one, 

M(6 2 ) = (A 2 | ri (A 3 | r2 ~ |ri ~ r21 |A 3 ) ri |A 4 ) r2 
can be derived from (|57]). Thus we obtain 
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and arrive at eq. (102). With the aid of this formula we find 
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With the aid of eqs. fllQl]) and jl03|) we observe that the result® is zero at T -> and it is equal to (5cr^ anh jlO^ ) 
for vrTt > 1. Combining eqs. (g^ @, @, @, |L0l]) and (|103D we arrive at the final results for <kr tanh : 



5a 



tanh 



7T (71 



1 Dt 



3/2 



2£r H 3C(3/2) /DVjf 



iry/2 y/T e it 2 



4tt 



2ttT 



itTt h -C 1, 



(104) 



for 



tanh 



7T (7i 



1 L»T 



3/2 



Tr-y/2 



CO/2) 
2\72^ 



dVTt^ /2 - 



3C(3/2) /£ 2 r ff 



2tt 



2ttT 



ttTt h > 1. 



(105) 



In order to find the total expression for the weak localization correction one should simply add the two contributions 

3/2 



Sa dcph @, (|f) and <5cr tanh (|oJ), ( ^05|) together. We observe that the temperature independent terms 



are equal in these two expressions, they enter with the opposite signs and cancel each other exactly in the sum 
<5<7 deph + <5i7 tanh in both limits irTt <C 1 and irTt 3> 1 . As we have already discussed, these are just the linear in time 
"golden rule" terms coming from the exponent (<5cr deph ) and the pre-exponent (<5<7 tanh ). Their cancellation occurs in 
no relation to (and due to much more general reasons than) averaging over disorder. Other ("non-golden-rule") terms 
do not cancel and combine in the final result which we will present below. 



D. Discussion 

Although the main differences between our approach and that of AAgE! can already be understood from the above 
analysis, we will briefly summarize them again for the sake of clarity. 

1. The first crucial difference tO|-be emphasized here is that our methodi is essentially nonperturbative in the 
interaction while the approacfiij is only the first order perturbation theory. In the most interesting limiting case 
tjj ^ (which was only considered in our Refs. one cannot proceed perturbatively in the interaction at 
any temperature including T = O.i-This is precisely what A AG do: it is demonst rated in Appendix A that the 
general result for the conductivityli3 is identical to the first order expansion of (A16) in the interacting terms 



iSn + Si whiler-sJl higher order terms (which are larger than the first order term for t# > t v ) were not taken 
into account int2l. In contrast, our path integral approach is equivalent to an effective summation of diagrams 
in all orders with the exponential accuracy. This is sufficient for correct evaluation of t v . Within our analysis 
only the action in the exponent (rather than the whole expression for a) is expanded in the interaction (this 
is correct as long as pfI 3> 1). Our method also allows for a clear distinction between the exponent and the 
pre-exponential contribution to 5a. 

It remains -unclear to us why AAG repeatedly stated, that our procedure "is nothing but a perturbative 
expansion"Ed and our results are "purely perturbative"llj. The only justification of the above statements which 
we could extract from the above papers is that our result for the dephasing rate 1/t v "is proportional to the first 
power of the fluctuation propagator" Ej. Although the latter is true in some limits, it is hard to understand how 
this could help to turn a nonperturbative problem into a perturbative one. Indeed, if one formally multiplies the 
photon propagatoi-by a constant A everywhere in our calculation, one would obtain 1/t v oc A. The same holds 
for the calculation^. Note, however, that it is not the decoherence rate 1/t v (we are not aware of a quantum 
mechanical operator which expectation value would correspond to such a quantity) but rather the expectation 
value for the current operator which is calculated theoretically and measured in experiments. For th ^ t v the 
result for the weak localization correction 5a depends on A as 

5a\ cc — 1/vA 
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in Id (cf. eq. (]{])) and S<j2 oc In A in 2d. Obviously, these results are purely nonperturbative in the "interaction 
strength" A. Any attempt to calculate the expectation value of the current operator perturbatively may only 
yield to divergences in all orders of the expansion in powers of A. As to the decoherence rate 1/t v , it is only 
extracted from the nonperturbative results for the conductance correction. Hence, the relation 1/t v oc A cannot 
by itself tell anything about the perturbative or nonperturbative character of the calculation. In the limit 
t~h •C Tip the conductance correction 5a{H) can be evaluated perturbatively in A. However, as it was explained 
in Sec. 2, even in this limit t v can be unambiguously determined only within the nonperturbative procedure, 
while any perturbative expansion yields ambiguous results for t v which fully depend on the assumption about 
the decay of correlations in time. 

2. Another crucial difference is that AAG essentially use the assumption about a purely exponential decay of the 
phase correlations in time while no such assumption was used within our analysis. Specifically, eq. (3.2) of Ref. 
|l3| is equivalent to our eq. (|^) only provided one assumes that fd(t) is a linear function of time fd(t) — t/r v 
and ignores the effect of the interaction on the pre-exponent, i.e. puts Adit) — 1/ (A-n:Dt) d l 2 . This assumption 
cannot be checked within the perturbation theory in the interaction and, as it was already explained above, in 
general it can only be valid within the golden rule approximation. The whole comparison between ours and 
AAG's results carried out by the authorsfca is essentially based on their eq. (3.2) which was neither used nor 
even written down in our papercl. 

Let us emphasize that AAG (unlike many others) do not use the golden rule approximation in their perturbative 
calculation of the weak localization correction Sa(H) in the limit th <C t v . However, they explicitly use this 
approximation while extracting t v from 5a(H): eq. (3) of Ref. [ll] and eq. (4.3) of Ref. 13 are valid only within 



the golden rule approximation. As it was demonstrated above, fd is not a linear function of time (cf. eqs. 
(|33|)) and, moreover, in the presence of interaction the pre-exponent in (||) deviates from its "noninteracting" 

form Ajj (t) = l/(4irDt) d / 2 . As a result, the relation between Sa(H) and t v depends on temperature and is 
different from eq. (3) of Ref. [ll] (or eq. (4.3) of Ref. |l3|) at any T even in the limit th -C t v . Since the linear in 
time T-independent contributions from the exponent and tha,ptp-exponent exactly cancel each other in the first 
order perturbation theory, the golden-rule-type assumptiontlEJ about purely exponential decay of correlations 
in time inevitably yields to missing of the T-independent contribution (^2|) to t v . 

3. Let us compare all the approximations used by AAG0 and in our paperi. In both papers the same quasiclassical 
condition pfI 3> 1 was assumed and the expressions for the photon propagators were defined within RPA. In 
order to perform the perturbative expansion in the interaction AAG considered the limit of strong magnetic 
fields th <C t v (in Ref. [ll] this condition was not quoted). AAG also performed the expansion in the inverse 
dimensionless conductance l/g(Ln), i.e. they assumed that 5 > 1 on the scale of the magnetic length Lh = 



\/Dth- Although we do not rp<;d these approximations within our nonperturbative analysisa, their appearance 
in the perturbative treatmentliia is understandable. 

As to an additional condition Tth 3> 1, in our opinion it is not needed even within the perturbative procedure 
of AAG. Indeed, the condition g(Ln) 3> 1 does not depend on temperature at all, and the inequality th <C t v 
can only become stronger at lower T provided it is akeady satisfied at higher temperatures. Therefore under 
the two latter conditions the perturbative expansiontia should be justified down to T — > and the condition 
Tth S> 1 is not needed at all. This condition should also be irrelevant for eqs. (2.42) of Ref. [ll| According to 
AAG "all the corrections to these formulas are small as 1/(Tt v )". Combining th <C t v and g(Ln) 3> 1 with 
eq. (2) of Ref. [ll] (or eq. (4.9) of Ref. ||) 1/t aag = T/g{L H ) we observe that the inequality Tt aag > 1 is 
satisfied at all temperatures including T — > 0. 

Thus the perturbative resultsEl can be analyzed in both limits Tth > 1 and Tth < 1- Since the latter limit 
of lower tempera tures was n ot discussed by AAG we carried out the corresponding analysis in Appendix C. 
Combining eqs. (|C16|) , ( |C17|) with eq. (4.3b) of Ref. [l| we find 1/t aag = Te 2 ^/Dt~h~ / {Aox) for Tt h > 1 (just 



like in Ref. [D| and 1/t aag = 3e 2 v r Dr H ' / \2-kgith) for Tt h < 1. The latter result (which was not presented 
by AAG) demonstrates that a nonzero dephasing time at T — is obtained even if one explicitly follows the 
procedure of Ref. [l3[ Although due to the reasons explained above this result differs from the correct one (|3^ ) 
it is interesting to observe that a nonzero dephasing rate at T = is already contained in the formulas derived 
by AAG. 

4. Subtle details of disorder averaging do not play any significant role in the problem in question and can merely 
influence some numerical prefactors of order one. As it was demonstrated above without making any approxima- 
tion, no exact cancellation of the first order diagrams occurs even at T = 0. The "non-canceled" T-independent 
terms describe not only renormalization due to interaction but also contribute to dephasing. These conclusions 
are general and hold both before and after averaging. 
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As to the (partial) cancellation, it indeed occurs at T — > 0, but only for the linear in time "golden rule" 
terms coming both from the exponent and the pre-exponent. This partial cancellation is also due to very 
general reasons, it occurs already in the exact (non-averaged) perturbative expression and has no relation to 
the quasiclassical approximation and/or disorder averaging. AAGO argued that this cancellation can not be 
reproduced if averaging over disorder does not involve the so-called Hikami boxes. Our analysis demonstrates 
that this is not true. 

Let us compare the perturbative results for the weak localization correction obtained in Ref. [l^ and within our 
analysis. In the limit Tth 3> 1 for the Id case AAG get (cf. eq. (4) of Ref. [n] or eq. (4.13a) of Ref. [lg| ): 
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In the same limit with the aid of our eqs. 
Sa tanh we find 
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In the opposite limit Tth 1 our calculation of the integrals!^ (see Appendix C) yields 
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Combining our eqs. ( p2| ) and (104) in the same limit Tth < 1 we obtain 
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Note, that the "renormalization" terms oc WthJ T (wh ich are irrelevant for dephasing and can be added to the 
interaction correction) are dropped in eqs. (106)- (109) for the sake of simplicity 

We o bserve that in both limits Tth 1 and Tth <C 1 the T-independe nt te r ms oc t^ 2 (see eqs. (|S2"|), ([83"|), ( 104 ), 



(105)) exactly cancel each-Qt her a nd do not contribute to the results (107), (10£) at all. The same cancellation 
occurs in the expressionstU (10£), (108). The latte r equa tions were derived within the averaging procedure 
involving the Hikami boxes. In order to obtain (|l07| ), (109) we used a somewhat different averaging procedure 
which amounts to deriving the matrix elements M(lj') from the general properties of diffusive trajectories. Since 
in both cases exactly the same cancellation occurs in both limits of high and low temperatures, we conclude 
that the issue of the Hikami boxes raised in Ref. [l^ is completely unimportant for this cancellation. 

We can also add that the averaging procedure employed by AAG is efficient within the perturbation theory 
while our procedure is developed to average the nonperturbative results obtained within the path integral 
technique. The perturbative results obtained within both methods are essentially the same in 2d (see Sec. 3C) 
and practically the same in Id apart from some unimportant details. Both procedures yield nonzero "dephasing" 
terms even at T = 0. 



V. CALDEIRA-LEGGETT MODEL 

As we have already discussed befor JB@0, the physical nature of the interaction-induced decohcrence can ..he 
understood with the aid of a simple model of a quantum particle interacting with a bath of harmonic oscillatorsHEJ. 
By a proper choice of both the interaction term and the frequency spectrum of the bath oscillators one can easily 
realize the important limit of Ohmic dissipation and arrive at the Caldeira-Leggett (CL) modeB. Some rigorous 
results obtained within this exactly solvable model are presented in Appendix D for the sake of completeness. 

An important advantage of the CL model is that the density matrix and the expectation values of the quantum 
mechanical operators can be calculated exactly. This enables one not only to avoid worries concerning the validity 
range of various approximations, but also to test these approximations employed in some other models which cannot 
be solv ed , exactly. .. I ia-jparticular, here we are interested in checking the approximations which have led various 
authors^llJild'LaEZllla'Ea to the conclusion about the absence of interaction-induced decoherence in disordered metals 
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at T — 0, or to the conclusionllS that the quantum correction to the classical decoherence rate is small and decreases 
this rate below its classical value. Since it is well known that the off-diagonal elements of the particle density matrix 
p(x, x') are suppressed due to interaction with the CL bath even at T = in equilibrium (this effect is nothing 
but nonzero decoherence at T = 0), it is interesting to test if it is possible to reproduce this result within the 
approximations employed in the above papers. 

Also, it is sometimes speculated that the results derived within the CL model cannot be compared to ones obtained 
for electrons in a disordered metal because of different statistics. One could conjecture that electrons in a disordered 
metal should have zero decoherence rate at T — predominantly due to the Pauli principle, while in the CL model 
nonzero decoherence at T = is allowed because no exclusion principle exists for bosons. The role of the Pauli 
principle can also be clarified by performing a direct comparison of the results obtained within the CL model with 
ones for electrons in a disordered metal. 

On a perturbative level this program will be carried out in the subsection A. In the subsection B we will discuss the 
relation between the exponent and the pre-exponcnt for the CL model and illustrate the analogy between the results 
of this subsection and those of Sec. 3. We will develop this comparison further in the subsection C where we analyze 
the properties of the "Cooperon" in the CL model. The validity range of various approximations is discussed in the 
subsection D. 



A. Perturbation theory 

Since in practically all cases the conclusion about the zero decoherence in the interacting systems at T = in 
equilibrium was reached only within the first order perturbation theory in the interaction, it is instructive to examine 
the structure of the first order perturbative terms in the CL model. 



Let us expand the kernel of the evolution operator (Dl) in the interaction part of the action iSn + Si. In the zeroth 
order we get a simple result 

Jo(t,Xlf,X 2 f,X ll ,X2i) = U(t,Xlf,Xu)U + (t,X 2 f 1 X 2 i), (110) 

where U(t, Xf, Xi) — (xf\ exp(— ip 2 t/2m)\xi) is a free particle evolution operator. Investigating the transport properties 
of disordered conductors one usually expresses the results in terms of advanced and reta rded Green functions G R ' . 
In order to emphasize the analogy with the CL model, we note that the expression ( |ll0|) can be rewritten as 

J (t, Xif,X 2 f 1 Xu,X 2 i) = G R (t, Xif,Xu)G A {-t, X 2 f,X 2i ), (111) 

where G R (t,xj, Xi) — —i9(t)U(t,Xf ,Xj), G A (t,Xf,Xi) = i0(—t)U(t,Xf,Xi). Comparing this expression to that for the 



conductivity of a disordered metal (A15), we note that the latter contains an additional time integral, a oc J dt Jo(t). 
This difference is not important though, in order to simplify the comparison of the corresponding perturbative results 
one can always keep the time t finite (exactly as it was done in the preceding section) and perform the time integration 
only at the last stage of the calculation. 

Let us consider the first order correction to the kernel J due to the interaction. This correction is again given 
by the sum of the four diagrams of Fig. 2. The current operators j a are, however, not applied. Also the "photon 
propagators" are now different. Namely, instead of the function i?(<3 — £4, — r±) one should subst itut e the function 



®R(t3 — t^x^x^, while instead of I{t^ — £4, r% — r^) one should use oti(tj, — t^x^x^ (see eqs. (D4), (D5)). In contrast 
to the case of an electron propagating in a disordered metal ( Al6- |Al9f ) the action in the exponent of (Dl) does not 



contain the factor 1 — 2n(j>. r). Therefore the operator 1 — 2p, related to the Fermi statistics, does not appear in the 
perturbation theory. The free particle states are labeled by its momentum, therefore the indices Xj in the diagrams 
of Fig. 2 should be understood as the momentum values. 

For the sake of brevity we will omit the general result for the first order correction S to the operator J which 
is expressed in terms of the same functions ( |58| ) and ( |6l| ) with £\ 3 — * E p . Rather we immediately go over to the 
part of the kernel J describing the evolution of the diagonal elements of the density matrix, which corresponds to the 
"diagonal" part of the conductivity (|64[). This is sufficient for our illustration purposes. For the probability of the 
transition from the state with the momentum q to the state with the momentum p after the time t we find 

SJ pl]qq( t ) = (pU lf (q\x 2l SJ {V ' ] (t,Xlf ,X U ,X 2 f , X 2 i)\p) x 2f \q) Xli = 

„ r 1 19 f du> ( 1 w \ 1 - cos \(E a — Ek — uj)t] 
-2 V S pq £ M » J _ W (coth - + l) {E *_< Ek _ u)2 

l9 f du ( . ui \ 1 — cos \(E V — E„ — , . 

+M*»?j ^(coth--l) {E *_ P Eq _l )2 } \ (112) 
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where x p k = (p\x\k) is the matrix element of the operator x and E p = p 2 /2m > is the energy of the free particle 
with the momentum p. 

Let us first evaluate the above expression within the standard golden rule approximation performed e.g. in Refs. 
|26|,|l7],|l8] and others. As it was already discussed in Section 4.2, this approximation is equivalent to the following 
replacement 



1 - cos [(E p — E q — Uj)t] 
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Performing this replacement in (|112| ) we get 
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(114) 



The first term in this expression, as well as in the expression (112), originates from the self-energy diagrams (a) and 
(b) in Fig. 2. This term describes the out-scattering rate. It is equal to the sum over all possible transitions from 
the initi al state q to all finite states k. It is mostly important as far as dephasing is concerned. The second term in 
the cq. (112,114) just gives the transition rate from a given initial state q to a given final state p. It comes from the 
vertex diagrams (c) and (d) in Fig. 2. 

Now let us investigate the time evolution of the density matrix provided the initial density matrix is just the 
equilibrium one for a free particle. In the momentum representation it has the form: 
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where L is the length of the system. Substituting all these results into eq. ([]]) we obtain the following expression for 
the occupation probability of the state p within the golden rule approximation: 
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We observe that the combination 



coth 



E p - E k 
2T 



(116) 



appeared in the result (115). It is very similar to the combination (|62j) in the case of a disordered metal. Keeping 
in mind the condition E p > 0, one can easily see that this combination again yields zero result at T = 0, i.e. within 
the golden rule approximation relaxation processes are forbidden in the zero temperature limit. This result is fully 
equivalent to one presented e.g. in eqs. (1,2) of Ref. and it is in an obvious contradiction to the exact solution of 
the CL model. 

Now let us perturbatively find the occupation probabilities without making the golden rule approximation. We get 
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As it was already done for the case of a disordered metal (Sec. 4B), let us consider the first part of this expression 
determined by the self-energy diagrams of Fig. 2a,b. In the zero temperature limit we find 
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where we defined E p k — E p — Ek- We observe a close similarity between the results (118) and (|67|). In both cases 
the first two terms in the right hand side are the same as in the golden rule approximation, they cancel each other 
at T = 0. In both cases the third term survives even at T = 0, it is due to quantum noise and originates o nly from 
the coth-part of the effective action. Moreover, we observe that the last terms in the equ ation s ( p7\ ) and (118) are 
exactly the same, one should only identify the energy difference £31 in ( |67| ) with E p k in (US). Thus we conclude 
that at T = the only difference between the two problems considered here lies in the matrix elements of the 
interaction. Everything else is the same and, hence, wc explicitly demonstrated that the Pauli principle cannot cause 
any important distinction between the problems in question at T = 0. The above difference in the matrix elements 
results only in some quantitatively different features, like e.g. different functional dependences of the density matrix 
on time, however the decay of the off-diagonal elements of p (and thus decoherence) is present in both models at any 
temperature including T = 0. In both cases at low temperatures this decay cannot be correctly described within the 
golden rule approximation. This approximation fails completely at T — > 0. Dropping the cos-term while evaluating 
the result (117) (this approximation is be equivalent to one we discussed in the end of Sec. 4B) is clearly insufficient 
at low temperatures (see below). 

Finally, a close similarity between the perturbative results obtained here and in Sec. 4B demonstrates again that 
averaging over disorder is absolutely irrelevant for the issue of cancellation (or non-cancellation) of the diagrams in the 
first order of the perturbation theory. In the CL model no such average exists at all, however at T = diagrams cancel 
or do not cancel depending on whether to employ or not to employ the golden rule approximation. This property is 
completely general, and it does not depend on the form of the matrix elements, Fermi or Bose statistics, averaging 
over disorder and other details of the model. 



B. Exponent and pre-exponent 



As the structure of the first order result in the perturbation theory is clear from the above consideration one can 
try to proceed further and calculate higher order terms. Then one can try to sum up diagrams in all orders in order 
to recover the nonperturbative result. This program appears to be quite involved from a technical point of view and, 
to the best of our knowledge, was not yet carried out. Fortunately in the case of the Caldeira-Leggett model one 
does not need to sum up diagrams, the exact result ca n be obta ined much easier, just by performing several Gaussian 
integrations. In this way one arrives at the equations (D6-D10) which are, of course, equivalent to the result of exact 
summation of all diagrams of the perturbation theory. The result (D6) can be expressed in the following form: 
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(119) 



where 7 = rj/m, the function R = R(t,Xi,Xt) does not depend on temper atu re , wh ile the function I(t) is propor- 
tional to the frequency integr al of t he co mbination u> coth(w/2T) (see eqs. (D6)-( D~10| )) rather than the combination 
w(coth(w/2T) + 1) ( see eqs. (pT^)-(|il7|)). 

By looking at the above formula (or at ( |D6D ) one immediately observes that this result cannot be simply guessed 
from the first order perturbation theory, e.g. just by exponentiating the first order results or by a similar procedure. 
As it was demonstrated above, the combination "coth + 1 " (|116| ) appears in the perturbation theory, while the time 
dependence of a real part I(t) of the exponent ( |l!9| ) (or (D6)) is governed only by "coth" and not by "coth + 1". 



This implies that in the course of the exact summation of all the diagrams terms combine in a nontrivial way, so 
that "coth" gets split from "1" , this combination does not appear in higher orders in the same form as in the first 
order result. From the exact result (D6)-(D1C) one can immediately draw a conclusion on how "coth" gets split from 
"1": all terms of the diagrammatic expansion cont aining "coth" gather in the exponent. Moreover, only such terms 
contribute to the real part of the exponent I ( |ll9|) . In other words, since "coth" is contained only in the imaginary 
part of the effective action Si (D3), we conclude that the real part of the exponent / (119) (or ( |D6| ) is determined 
only by Si. This result could, of course, be expected in adva nce because it is the imaginary part of the action Si 



which should be responsible for the decay of the kernel J (D6) in time. 

The real part of the action Sr (D2), in contrast, does not (and cannot) contribute to the real part of the exponent 
/. Now we will demonstrate that Sr (and only S r) c ontributes to the pre-exponential function in (119), (D6). The 
pre-exponent is determined by th e path integral (Dl) with zero boundary conditions, xn — x-n = Xif = X2f = 0. 
From the structure of the action (Dl-D3|) one easily observes that its real part can be represented in the form 



iSolx^ - iS a [x 2 } - iSR\x 1 ,x 2 ] = ix Lx + , 



(120) 



while its imaginary part has a different structure: 
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Si[xi,X2\ = x Ax , 



(121) 



where L and A are two different operators. The Gaussian integrals can be easily performed and we get 

f Vx- [ x>x + e lx ~ Lx+ - x ~ Ax ~ = [ Vx-6(L T x-)e- x ~ Ax ~ = — = ; V --. (122) 

J J J y ' detL 27r(l-e-T*) v ' 

This relation proves that the pre-exponent in ( |119| ) (or (p6|)) is determined solely by the real part of the action Sr 
and does not depend on its imaginary part Si at all. 

Thus the above analysis allows for a clear distinction between the two parts of the effective action: Si determines 
the real part of the exponent I(t) and governs the decay of the off-diagonal elements of the density matrix (hence, 
playing a crucial role for dephasing), while Sr determines the time dependence of the pre-exponent which is only 
relevant for the kinematics of classical trajectories and has nothing to do with the issue of decoherence. In the first 
order perturbation theory the terms from the exponent and the pre-exponent mix and partially cancel each other, 
thus making any clear distinction between them impossible. This situation is fully analogous to one encountered in 
the preceding sections for the problem of electron dephasing in disordered conductors. 

Depending on the boundary conditions Sr can also contribute to the imaginary part of the exponent R (cf. ( p6] ). 
However, this imaginary part does not determine the decay of the off-diagonal elements and in no way can it cancel 
(or contribute to) the real part of the exponent I(t) determined exclusively by Si. On top of that, for the boundary 
conditions corresponding to the "Cooperon" this imaginary part of the exponent vanish R = 0. This particular case 
will be considered below for illustration. 



C. "Cooperon" in the Caldeira-Leggett model 

Let us define the "Cooperon" configuration as 

C(t,x) = J(t, x, 0; 0, x) 



(123) 



That means the path integral (Dl) is evaluated on the trajectories with the time reversed boundary conditions: 



xu = %2f — 0, xif = X2i — x or, equivalently, x\ = Xj = x/2, x 
conditions into (|D(]) one easily finds 
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This is the exact result. In the long time limit jt 3> 1 we obtain 

C ~ 2- e -nai (*)* 2 ) gi ( t ) ~ Tt + ln 
2ir 



-x. Substituting these boundary 
exp [-v{9i(t) + 4s 2 (t) - 2g 3 (t)}x 2 ] (124) 
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This expression decays in time, and the analog of the dephasing time (which we will also denote as t v here) can be 
defined from the equation 



V9i(t v )x 2 w 1. 



(126) 



We observe that this time is determined only by the exponent of the equation (124). 

Let us no w expand separately the combination of ^-functions in the exponent and the pre-exponent in the exact 
result (124) to the first order in the interaction parameter ij each. We get 
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(128) 



The function g pcr t is just the imaginary part of the action Si(t) (D3) evaluated on the two time reversed paths: 
xi(s) — xs/t and ^(s') = x(t — s')/t. These are the saddle point paths for the noninteracting part of the action Sq. 
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Comparing the functions <? pC rt (128) and g%(t) 
for Tt < 1 and g pcrt (t) ~ gi(i)/3 for Ti > 1 



result (127) as 



123 ) we observe that with a sufficient accuracy one has g P ert(t) — gi{t) 
Hence, if the dephasing time t v is determined from the perturbative 
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(129) 



the result will differ from the exact one only by a numerical factor of order one which anyway can be absorbed in 
the definition of t v . As before (cf. ( |126| )) the dephasin g tim e extracted from this equation will depend only on the 
exponent (i.e. on Si), while the pre-exponcnt in the eq.(127) (defined by Sr) must be ignored again. 

Now let us expand the whole expression (124) in powers of -q. The first order correction for the Cooperon has the 
form 



8d(t,x) = C (t,x) 
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Wpert(i)x' 



(130) 



where Co(t) = rajl-nt is the Cooperon in the absence of the interaction defined e.g. by (111). The equation ( Jl 30| ) is 
nothing but the short time expansion of the exact result because the interaction parameter 77 always enters into this 
result being multiplied by t. This expansion is fully equivalent to the perturbation theory in the interaction in the 
weak localization theory. 

As it was already discussed in Sec. 2, it is impossible to unambiguously define the dephasing time t v within the 
frames of the perturbation theory only, and an additional assumption about the decay of the Cooperon in time should 
necessarily be made. For instance, AAG assumed that the Cooperon decay is purely exponential: C{t) = Co(i)e~ t//T<p 
(cf. e.g. eqs. (2.45) and (3.2) of Ref. [l3|). Assuming such a form here and combining it with the perturbative result 
(|l30|) , we find 
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(131) 



Comparing this equation with eq. ( |l2q ) which follows from the exact result (124) we observe the presence of an 
additional term —t]/2t: in eq. (131). This term originates from the expansion of the pre-exponent and has nothing 
to do with dephasing. However, it is not small and can strongly influence the result for t v , provided the latter is 
determined from the perturbative expansion ( |13C| ). Depending on the choice of x one can obtain positive, almost 
zero and even negative values of the dephasing time (cf. Sec. 2), which is an obvious nonsense. This simple example 
demonstrates again that it is impossible to make any conclusion about the long time behavior of the system (and, 
hence, about the dephasing time t v ) from the first order expansion in the interaction, since the latter is valid in the 
short time limit only. 



D. Other approximations 



Let us check some other approximations which are sometimes employed in the literature. A deficiency of the golden 
rule approximation applied within the perturbation theory has already been illustrated above. In certain situations 
one proceeds beyond the perturbation theory, correctly get "coth" in the exponent, and only then apply the golden 
rule approximation. Also in this case the true low te mper ature behavior will be missing complet ely. In order to 
observe this property let us eva luate the function gi(t) (D7) within the golden rule approximation ( 113 ). Extending 
the integral over s — s' in (D7) to infinite limits and pe rform ing this integration first we obtain the delta- function 
S(uj). This is just the golden rule approximation, cf. eq. ( |lI3| ). After that the remaining integrals trivially yield 

9l {t)~Tt, 



i.e. only the first term in (D7) (or (125)) is reproduced, while the second term is missing. After that one could 
incorrectly conclude that no quantum decoherence occurs in the CL model at T = 0. An obvious mistake here is to 
extend the integral over s — s' to infinite limits. An exact calculation of the function gi(t) allows to recover both the 
zero- frequency contribution Tt as well as an additional term (see eq. (125)) which originates from frequencies u> > T 
and does not vanish at T = 0. r— ■ 

The above approximation was applied e.g. in a recent preprint by LevinsonE3. In this paper a transparent formu- 
lation of the problem of quantum decoherence in quantum dots is derived. Within this formulation Levinson arrived 
at the result for the equilibrium dephasing rate which contains only "coth" , while "tanh" or "1" do not appear in 
the exponent at all (cf. eqs. (4), (5), (7) and (14) of Ref. [32]). It is interesting that a nonzero dephasing rate at 
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T = is contained in these formulas before the golden rule approximation is made. In order to see that it is sufficient 
to combine eqs. (5), (7) and (14) of Ref. [52], substitute them into eq. (4) of the same paper and perform the time 
integration. A finite result for the decoherence rate will follow immediately at all temperatures including T = 0. High 
frequencies will contribute to this result which will actually depend on a physical high frequency cutoff. However, if 
one applies the golden rule approximation, only the zero frequency contribution will be recovered (eq. (6) ofJicf. |3^ ) 
and the whole decoherence effect at low T will be missing. At low T the decay of correlations in the problen£3 is not 
exponential in time. But also in this situation the dephasing time t v can be easily defined. This time just sets a scale 
on which the quantum coherence is sufficiently suppressed and the integral over time (eq. (4) of Ref. ^) becomes 
convergent. 

For nonzero (but possibly relatively small) values of T one can try to argue that at long times it is sufficient to 



consider the limit Tt 3> 1 and expand the function g%(t) (D7) (or (|125|)) in powers of 1/Tt. Then one gets 



gx(t) « Tt + 1n(u c /T), (132) 
while all higher order terms of this expansion will be equal to zero because gi(t) is a non-analytic function of 1/Tt. 



From (132) one could conclude that the first term in this equation describes dephasing due to thermal fluctuations 
while the second term is the "interaction correction" which does not depend on time and has nothing to do with 
dephasing. Since the first term vanishes in the limit T — > 0, one could again arrive at the conclusion that no 



decoherence is present in the CL model at T — > 0. As it is clear from the exact solution (D6-D10), this conclusion is 
not correct. 

The expansion in 1/Tt (or in 1/Ttjj) is just the expansion performed in Ref. [l] within the first order perturbation 
theory in the interaction (see also Sec. 4D and Appendix C). As it was already discussed in Sec. 4B, this exp ansio n is 
equivalent to dropping the oscillating cos-term e.g. in the expression ([TO]) (or, equivalently, in the expression flllTj ) for 
the CL model). Obviously this approximation has nothing to do with averaging over disorder. The above example also 
invalidates any attempt to approach the correct low temperature behavior by means of a high temperature expansion, 
like it was suggested e.g. by Vavilov and Ambegaokar in Ref. [l(| For instancftHthe terms ~ exp(— 2itTt) cannot be 
recovered in any order in 1/Tt. We will come back to the analysis of the paperE3 in Appendix E. 

Finally, let us briefly discuss the role of the low frequency modes of the effective environment. One could conjecture 
that, since high frequencies up to lo c contribute to the dephasing time, the result will not change even if one introduces 
the low frequency cutoff ui c0 . In order to test this conjecture it is again sufficient to consider the behavior of the function 



gi(t). Introducing the low frequency cutoff u c o in the integral (D7), at T = one readily finds 



gi(t) ~ ln^c/wco), w c0 t>l, (133) 

while in the opposite limit Lo c ot < 1 the result is the same as without the low frequency cutoff, i.e. gi(t) ~ hiuj c t. 
Thus the cutoff at u> ~ ui c o leads to a different long time behavior of gi(t). It increases at short times but then 
saturates at a value ~ ln(w c /w c o). In this case no time decay of the off-diagonal elements of the density matrix occurs 
at long times and therefore the coherence is not fully suppressed. However, it is not the long time limit which is 
i nter esting in the dephasing probl em, but rather the system behavior at t ~ t v . If uj c qt v <C 1 (t v was defined e.g. in 



(126)), by the time the behavior (133) is reached the coherence will already be very strongly suppressed. Thus from 
a practical point of view there is no substantial difference between the cases lo c q = and lo c qt v <C 1. In the opposite 
limit u c oT<p 3> 1 the function gi(t) saturates earlier than quantum coherence gets suppressed. In this case cutting out 
the low frequency oscillators changes the result significantly. 

This simple consideration clarifies the role of the low frequency modes in the dephasing problem. At low T these 
modes do not really affect the expression for t v which depends on the high frequency cutoff uj c . However, if the 
low frequency cutoff is chosen such that lo c qt v 1, the dephasing time t v simply looses its meaning because gi(t) 
saturates already at much shorter times t ~ l/u> c o ^ T <p- The same conclusion applies to disordered metals in which 
case the function fd(t) (p3) should be considered instead of gi(t). 

For an extended discussion of various approximations analyzed for the exactly solvable CL model we refer the reader 
to the paperEa. 



VI. CONCLUDING REMARKS 



Let us briefly summarize our main observations. 

1. We explicitly demonstrated that the perturbation theory in the interaction is principally insufficient in the 
problem of electron dephasing in disordered conductors. The decoherence time t v cannot be extracted even 
from a correct perturbative calculation unless one assumes some particular form of the decay of correlations in 
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time. However, this form should be calculated rather than assumed. This task can be accomplished only by 
means of a nonperturbative calculation. 

2. A nonperturbative analysis shows that the dephasing time t v is determined only by the function fd(t) in eq. 

while the effect of interaction on the pre-exponential function Ad(t) is not important for the calculation of 
t v . Therefore, in order to evaluate the dephasing time it is sufficient to perform a nonperturbative analysis with 
the exponential accuracy. 

3. The effect of interaction on the pre-exponent Ad(t) is important if one calculates the interaction-dependent 
part of the weak localization correction in the limit of strong magnetic fields th <C t v . The zero temperature 
dephasing time r° drops out of this correction in the first order due to the exact cancellation of the linear in 
time T-independent contributions from the exponent exp(— fd(t)) and the pre-exponent Ad(t). 

4. Nonlinear in time T-independent contributions do not cancel out already in the first order of the perturbation 
theory. In general these terms not only account for the renormalization effects but also contribute to dephasing 
at all temperatures including T = 0. 

5. We demonstrated that there exists a close formal similarity between the problem of electron dephasing in 
disordered conductors and the exactly solvable Caldeira-Leggett model for a particle interacting with a bath of 
harmonic oscillators. 

Our analysis allows to establish a simple correspondence between the results of Refs. ||J| and |l3|. The effect of the 
interaction in the expression for the magnetoconductance ([}]) is described by the function 

A d [tanh] exp(-/ d [coth]). 

At high temperatures the function fd decays on a short time scale. For such times the effect of the interaction (related 
to "tanh") on the pre-exponent is negligible. In this case the nonperturbative analysis of AAKd applies. At low 
temperatures the interaction effect on both the exponent and the pre-exponent becomes of order one on the same 
time scale ~ r°. Since the effective cutoff in the integral (|^) is determined by the function fd(t) in the exponent, "tanh" 
can be neglected again. We arrive at our nonperturbative results!! Finally, in the limit th <§C t v for the relevant 
times t < tjj one has fd(t) <§; 1 and, hence, exp(— /^[coth]) ~ I — /^[coth]. Performing also a short time expansion-af 
yl^tanh] we obtain the combination "coth - tanh" in the first order and reproduce the AAG's perturbative resultsliil 
These observations conclude our analysis. 

We acknowledge stimulating discussions with B.L. Altshuler, C. Bruder, M. Buttiker, K. Eriksen, P. Hedegard, D. 
Khveshenko, V.E. Kravtsov, D. Loss, Yu. Makhlin, A. Mirlin, D. Polyakov, G. Schon, S. Sharov, E. Sukhorukov and 
J. von Delft. One of us (D.S.G.) also acknowledges the support from the INTAS-RFBR Grant N 95-1305 and the 
hospitality of the Karlsruhe University and the Forschungszentrum Karlsruhe where a part of this work was performed. 



APPENDIX A: 



Here we will compare the expressions for the weak localization correction to the conductivity in the presence of 
interaction obtained by means of our path integral techniquea and within the diagrammatic approach of AAG. We 
will demonstrate that both results are identical if analyzed on a perturbative level before disorder averaging. 

We proceed in two steps. We first transform the resulted for the conductance and demonstrate that by virtue of the 
causality pri ncipl e one can completely remove the terms of the tyn&jG R G A G R G A emphasized by AAG. We will arrive 
at the Eqs. (|All| , JaI| ) wh ich are exactly equivalent to the resultt^l. Our second step is to ex pand our expression for 
the cond uctance (A15,A16) in the interaction terms iSa + Si. This will lead us to the Eq. (A23) which is identical 
to© 

We start from reproducing the expression for the correction to the conductivity due to electron-electron interaction 
obtained in Ref. [L3] 
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{2j a [Gf 2 (e) - G A 2 (e)] h [G&(e)Gft(e)r&M - G«(e)G« (e)£ 3 R 4 (o,)] [G* (e - w) - G&(e - u,)] 

+... + « (Al) 

For simplicity we keep the same notations as in Ref. |lj: G R ^ are the retarded (advanced) Green functions for 
noninteracting electrons and C R ^ A ' are photon propagators. The coordinate dependence of the propagators is indicated 
by the subscripts, e.g. G R 2 {e) — G R (e, r\, r 2 ). Note that in (|Al| ) only the contribution of the two self-energy diagrams 
(Fig. 2a, b) was reproduced, while the remaining contribution from the vertex diagrams of Fig. 2c, d which contains 
the terms with two w-dependent Green functions is denoted by ... . We will not consider them in this Appendix for 
the sake of simplicity. 

We observe that the factor tanh enters in this expression together with the difference [G R A (e — uu) — G A 4 (e — w)l . 
This combination is just the Keldysh function 



G K (e, ri.rj) = tanh ^ [G R (e, n,r 2 ) - G A (e, n,r,)] = tanh ^ 



1 



This function can also be rewritten as follows 
G K (e,n,r 2 ) =tanh-i;^ 



e - f a + *0 e - £ A - *0 
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(A2) 



X)(-27ri) (tanh fA)*A(ri)^(r a ) 

E| tanh 
V 2T 
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#A(ri)* A (r 2 ) 



<*r' [G fl ( e ,ri,r') - G^e.n/)] (J(r ; - r 2 ) - 2p(r',r 2 )), 



(A3) 



where ^ are respectively the eigenvalues and the eigenfunctions of the Hamiltonian H — fi; p(r',r 2 ) is the 
equilibrium single electron density matrix, p = l/(exp((iJ — p)/T) + 1). In a similar manner one obtains 



tanh^) [G R (e,r l3 r 2 ) ~ G^(e,r 1 ,r 2 )] = 2 j dr' ^f^ [G R (e,r',r 2 ) - G A (e, r'. r 2 i 



(A4) 



We also introduce the evolution operator U(t) = exp(— i(H — p)t) which is defined both for positive and negative 
times. The functions G R and G A are related to this operator by means of the following equations: 



G R (t, n, r 2 ) = -iO(t)U(t, n, r a ); G A (i, n, r 2 ) = iO(-t)U(t, n, r 2 ). 



(A5) 



Now let u s wr ite down the two equivalent forms of the Keldysh Green function in the real time representation. We 
find from Q : 



+ 00 
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nT(t-f)) U{t '> ri ' r2) > 



(A6) 



and from ( |A^ ) we get 

G K (t, ri ,r 2 ) = J dr> [G R (t,r 1 ,r')-G A (t,r 1 ,r')] (8(r' - r 2 ) - 2p(r',r a )) 
= -*/ dr / t7(t,ri,r / )(«(r / -r 2 )-2p(r / ,r 2 )). 

Analogously we obtain 



(A7) 
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^tanh^) [G R (e,r u r 2 )-G A (e,r u r 2 )] => j <* ^ { P_ V)) [C*(f,n,r a) - G\i> >x,r a 



— OO 
+00 
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sinh(7rT(i - i')) 



(A8) 
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U(t,r',r 2 ). 



(A9) 



It is easy to observe that the eqs. (A6|,A^) contain the integral over time which does not enter the eqs. dA7L|A9|). It 
is this additional time integration that leads to violation of the normal time ordering at the level of the perturbation 
theory and is responsible for the appearance of the diagrams G R G A G R G A . The interpretation of such dia grams 
in terms of the path integral is not possible. However, if one uses the other form of the same functions ( A7, A9j) 
the normal time ordering is automatically restored, the combinations G R G A G R G A disappear due to the causality 
principle and the path integral interpretation of the remaining terms of the perturbation theory can be made. 

We emphasize that all the above transformations are exact and have the advantage that in the final expressions 
only the propagators depend on the frequencies e and u> (except for the factor coth ^ in 6cr a p). This allows one to 
use the analytical properties of the propagators related to the causality principle. Namely, G R (e) and £ R (cu) have no 
singularities in the upper half-plane, while G (e) and £ (<jS) are analytic functions in the lower half-plane. Making 
use of these properties one can easily prove the identities 



(A10) 



du£ R (u)G A (e -w) = 0,J d e Gf 2 ( e )G 2 4 3 (e)G^ 4 (e - w)G&(e) = 0, 
J du£ A {u)G R {e -u) = 0,J deG R (e)G R (e)G R (e - cu)G R (e) = 0. 



Consider e.g. the integral J duj£ R (uj)G A (e — ui). Since both functions £ R (co) and G A (e — ui) are regular in the upper 
half-plane, the integral vanishes. Alternatively, we can write J du£ R (uj)G A (e — ui) = J dtexp(iet)£ R (t)G A (t) and 
note that £ R (t) = for t < due to the causality principle, whil e G A (i) = for t > and the integral is identically 
equal to zero. Analogously one can prove all the other identities (|A10 ). 
The corrections to the conductivity can now be considerably simplified: 
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We observe that the terms of the type j a G R 2 (e)j i3G 2 i 3 (e)G 34 (e — ui)G A x {c) do not enter the expression ( All ) at all. 
For later purposes it will be useful to rewrite the above expression in the form of the time integral: 
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where 



R ^ - 1 wFm'^ = ^ -h*™> (A13) 
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Now we will demonstrate that the equation (A.12) can be obtained within the path integral formalism. The formal 
expression for the conductivity has the fori 
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The kernel J is given by the path integral over electron coordinates and momentums ri(t),pi(t) and r 2 (t),p 2 (t) 
corresponding respectively to the forward and backward parts of the Keldysh contour. The explicit expression for 
this kernel readsQ: 

r 1 (t)=r lf r 2 (t)=r 2 f 



J(t,t';r 1 f,r2f,rxi,r 2 i)= J T>r 1 J Vr 2 J Vp 1 J Vp 2 x 



ri(t')=rij r 2 (t')=r 2i 

xexp{iS [ri,pi] -iS [r 2 ,p 2 ] - iS R [r 1 ,p 1 ,r 2 ,p 2 ] -5/[ri,r 2 ]}; (A16) 

where 

S [r,p}= f dt" Lr-E-- U(r)) ; (A17) 
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v v 
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+R(h - t 2 ,ri(t 1 ) - r 2 (i 2 ))[l - 2n(p 2 (t 2 ),r 2 (* 2 ))j 

-ii(ti - «2,r a (ti) - ri(t 2 ))[l - 2n(pi(t 2 ),n(t 2 ))]}; (A18) 
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The functions R(t,r) and I{t,r) are denned i n (|A13 ,A14) 



In order to obta in the perturbative result ( A12 ) from the formally exact expression ( A15 ) one needs to expand 
the kernel J ( A16| ) in iSr + Si. In the first order one obtains eight different terms. Again we will consider only the 
terms originating from the self-energy diagrams of Fig. 2a,b, i.e. the terms containing R(t\ — t 2 ,ri(ti) ~ fi(t 2 )), 
R(ti — h, r 2 (t\) — r 2 (t 2 )), I(ti — h, r i(ti) — r i{t2)) and I(t\ — t 2 ,r 2 (ti) — r 2 (t 2 ).) Four other terms which relate two 
different branches of the Keldysh contour and contain both r\ and r 2 , come from the vertex diagrams of Fig. 2c, d. 
As it was already discussed before these terms determine only a part of the function Sfd(t) ( |3l| ) and - although they 
do not vanish even at T = - always yield only a subleading contribution to fd(t)- Therefore we will not consider 
these terms here for the sake of simplicity. 

The correction to the kernel J due to the term I(ti — t 2 ,r%(ti) — ri(t 2 )) has the form 
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(A20) 



Here we made use of a simple property of a path integral: 

r{t)=r f 



J VrJ Vpf(t",r(t"))exp{iSolr,p]} = J dr"U(t - t";r f ,r")f(t" ,r")U(t" - t;r" ,r,), (A21) 
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which holds for an arbitrary function f(t",r(t")). Actually in deriving ( A20 ) the property ( A21 ) was used twice 
because the function of two arguments I(t 2 — t 3 ,ri(t 2 ) — ri (t 3 )) enters under the integral ( A20| ). Alrea dy at this 
stage one can observe the similarity between the expression (A20) and the second term in the expression (A12). To 
estab lish the equivalence between these two exp ressi ons the following steps are in order: i) after substituting the result 
( |A2C| ) into the expression for the conductivity ( ]A15[ ) and applying the current operator j — (ie/m)(V ri/ — V r2/ ) one 
puts rif — r 2 f = r 2 , ru — n, r 2 i — r$; ii) one denotes t — t' — > t\, t — t 2 — » t 2 , t — 1 3 — > ^3; hi) one introduces an 
additional integration J dr 2 /V which is just averaging of the expression (A15) over the sample volume and iv) one 
transforms the effective initial density matrix as follows 
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After jthcfc transformations one can immediately observe the equivalence of the results obtained by means of two 
methodsolij on the level of the perturbation theory. The terms arising from the real part of the action Sr can be 
transformed analogously, the only differ ence in this case is the presence of an additional factor (1 — 2,0)34 related to 
the term 1 — 2n(p, r) in the expression (AlS). Finally we get 
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(A23) 



In order to verify complete equ ivale nce of ( A12 ) and ( A2S ) one should a) replace the operator ep/m by j; b) adjust 
the factor 3 by observing that (A23) and (A12) are the corrections respectively to to the scalar conductivity and the 
conductivity tensor (in the i sotro pic case one has 5a = (Sa xx + Sa yy + 8a zz )/S) and c) adjust another factor 2 having 
in mind symmetrisation of ( A12Q with respect to indices a and (3. Also, one should keep in mind that the ope rator 



dpis/ d/i commutes with the evolution operator {7 52 . This completes the proof of equivalence of the results (A12) and 

The difference between the diagrammatic representations of the A AG's perturbation theory and in our approach 
is illustrated in Fig. 3. We have demonstrated that the initial diagram (Fig. 3a) can be expressed in two equivalent 
ways: by means of the diagram of Fig. 3b or as a sum of the two diagrams of Fig. 3c. In the first case causality is 
explicitly maintained and an additional coordinate integration appeared together with the density matrix 1 — 2p. In 
the second case only the sum of the diagrams of Fig. 3c is meaningful. Causality is violated in each of this diagrams if 
one considers them separately. This is particularly clear from Fig. 3d which shows the classical paths corresponding 
to the second diagram of Fig. 3£c) (cf. Fig. 10c of Ref. |l3|) According to Ref. 13 it is this path configuration which 
was "mistreated" in our analysisB. In Fig. 3d we observe that electrons move backward in time between the moments 
t' and t\. Such paths cannot appear within our path integral formalism, they are forbidden by the causality principle. 
Hence, their "mistreatment" could not occur within our analysis either. 




l-2p 




-Kftrf 





FIG. 3. The first order diagrams which contain the Keldysh function G K and tanh%^. (a) Initial diagram, (b) G K is 
written in the form (|iv?|). (c) G K is written in the form ( 

contains Hikami boxes. Here we have defined K(t) = sin ^ I J Tt ^ . (d) The classical paths which correspond to the second diagram 
of Fig. 3(c). Such paths violate causality and, hence, do not appear in the path integral. 



A6). Two diagrams appear. The second diagram G R G A G R G A 
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APPENDIX B: 



In this Appendix we will present some details of the derivation of eqs. (|5q-|6l|) . Let us first consider the contribution 
of the self-energy diagrams of Fig. 2a, b and consider only the terms in ( |A12[ ) which contain the function I(t, r), or 
coth(w/2T). Substituting eqs. (|53]-p5[) into these terms of (A12) we find: 
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e 2 f dridr2dr^dr4 
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1 / d 
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OO t\ t2 
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tanh ^ ) j a ^ Xl (r-OVl (r 2 )e~^j^ 2 (r 2 )^ 2 (r 3 )e^ 



2 \d£ Xl 2T 

[Ifa - *3, r 3 - r 4 )V-A3(r-3)^ 3 (r4)e^3fe- t 3)] ? /, A4 ( T . 4 ) ? /,* 4 ( ri )e^4(*i-^) + a «_> + c . c . 
This complicated expression can be rewritten in a simpler form if we introduce the matrix elements 



jAiAs;A,A a ^ = j dridr ^ (r x )^ (n)J(t, r x - r 2 )^ 3 (r 2 )^A 2 (r a ) 
and @. Then we get 
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Ai ..A 4 n n n 



(ta-*s) 



000 



COs(-fA 1 ti+fA a t3 + CA,(t2-*3)+fA 4 (tl-*2)). 



(B3) 



Analogously, one can find the contribution of the remaining terms in (A12) containing the function R(t,r), or 
tanh(£ A /2T): 
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J dtij dt 2 J dti 
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tanh sin (-£ Al *i + £ A2 t 3 + £a 3 (*2 - i 3 ) + - * 2 )) . 



(B4) 



Now we rewrite the functions / a 2A 3 ;A 3 a 4 ^ 2 _ ^ ancl ^>a 2 a 3; a 3 a 4 ^ 2 _ ^ ag f u ows 

^A 2 A 3 ;A 3 A 4 



l (*2-*3) 



do; d 3 fc / -47r 
Im 



27r(27r) 3 \k 2 e(u,k) 



coth — (A a | e ikr |A 3 ) (A 3 | e~ ikr |A 4 ) cos(o,(i 2 - f 3 )), (B5) 



R^^(t 2 -t 3 )=l ^^ F ^)(A 2 |e^|A3)<A 3 |e-^|A4)e-^3). (B6) 

Now let us express the kernel R ( B6 ) in terms of the imaginary part of the inverse dielectric susceptibility. Firstly we 
write 1/e = 1 — (1 — 1/e). The function 1 — 1/e is regular in the upper half- plane of u) and it tends to zero if u) — ► 00. 
Therefore for t 2 > £3 one has 

7 2tt ^e(w,fc) / 

Since in the integral (B4) the time t 2 indeed exceeds i 3 , we can replace e~ 4 "(* 2 ~ t3 ) by e ~M*2-*3) _ e M*2-i 3 ) m ^h e 
integral (B6). Then we arrive at the following result: 

flA 2 A 3 ;A3A 4(t2 _ h) = d{h _ t3 _ Q) (Aa | pi (A3 | r2 1 ^3)^ ^ 



/ 7*2 



2tt 



(B7) 
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The correction to the conductivity due to the self-energy diagrams can now be written as 

+00 



X se x C e 
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(A 2 |e Ite ^ |A 3 > <A 3 |. 



-ikr 



|A 4 ) 



coth — — + tanh 

2T 2T 



F(ti, (B8) 



where 5cr^o ([59|) is the correction due to the non-screened Coulomb interaction and the function F is defined in eq. 
Now let us make use of the Drude approximation for the dielectric susceptibility of a disordered metal 

47TCT 



e(o>, k) = — : ,. w ,, - 

v ; -lu + Dk 2 

In this case it is convenient to introduce the matrix elements M A2A3;AaA4 as defined in eq. (f37|). Combining the above 
expressions we immediately arrive at the final result for the first order correction to the conductivity from self-energy 
diagrams of Fig. 2a,b @, @. 

The conductivity correction from the vertex diagrams of Fig. 2c, d is evaluated analogously. After a straightforward 
algebra we obtain 
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Introducing again the matrix elements (B7j) and combining (B9) and (BIO) we arrive at the result (p0[), (61) 



APPENDIX C: 



In this Appendix we will analyze the expressions for the weak localization correction to the conductance obtained 
by AAG in the limit th <C t v perturbatively in the interaction. This additional analysis is necessary because the 
final results for the most interesting low temperature limit Tth ^5 1 were not presented in Ref. |l3|. For the sake of 
definitencss we will consider only the Id case which is sufficient for our purposes. 

AAG split the total expression for the weak localization correction to the conductance into two terms: 



5a 



AAG 
WL 



5a 



AAG 
deph 



5a, 



AAG 
CWLi 



(CI) 



where - according to Ref. [l3|- the first term describes dephasing while the second term accounts for the renormalization 
due to interaction and has nothing to do with dephasing. Below we will demonstrate that such a separation of the 
weak localization correction into "dephasing" and "interaction" terms is not possible even within the perturbation 
theory employed by AAG. 



Let us consider the first term in (CI) which was defined in eqs. (4.5), (4.6) of Ref. |l3| and has the form 

2 1 / AT ^2 
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4Te 2 1/0 

) [[C(Q,0)} 2 C(Q + q,co)-\C(Q 1 iu)\ 2 C(Q + q,0) , (C2) 



where the Cooperon and the diffuson are defined as follows: 
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The integrals over Q and q in (C2) can be evaluated exactly and we find 
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Making use of this equation we get 
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In the limit Ttjj ^> 1 the expression (C5) practically coincides with the analogous expression obtained in Ref. |l3j 
(see eq. (4.11) of that paper) with the only difference in a numerical coefficient in front of the last term in the first 
equation ( |C5| ) (5/8 in our calculation and 1/2 in Ref. 13). In the opposite limit Tth <C 1 the expression (C2) was 
not evaluated by AAG at all. 

Now let us consider the second contribution (CI) which was denoted by AAG as a cross term of weak localization 
and interaction. According to Ref. |l3| this term has the following structure (see eq. (5.23) of that paper) 
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where the terms in the square brackets arc defined in cqs. (5.25) o: 



I 2 (u> 
I 3 (uj 
I 4 (uj 
I 5 (oj 



dQdq 1 
(2tt) 2 Dq 2 
dQdq q d 
(2tt) 2 Dq 2 ~dq~ 
dQdq 

Wf 

dQdq Q 



C 2 (0,Q) 



C(-UJ,Q + q) 



V(-u,q) 
C(-uj,Q + q)C(0,Q), 



[C(-w, Q) - C(0, Q)] V(-lu, q) + 2C(-w, Q + q)C(0, Q) 



(27T) 2 q 
dQdq 



V(-oj,q)C(~OJ,Q + q)C(0,Q), 
V(-oj,q)C(0,Q-q)C(-UJ,Q), 



(27T) 2 

dQdq 
(2tt) 2 



V(-0J,q)C 2 (-UJ,Q), 
Dq 2 V 3 (-uj,q)C(-w,Q). 



(C7) 



These integrals can also be evaluated exactly and we arrive at the following results: 
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Let us compare these exact expressions with those obtained in Ref. 13. Unfortunately AAG calculated the integrals 
approximately only in the limit of high frequencies loth ^> 1 . In this limit they found 
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In order to find a numerical coefficient in front of the term ~ l/(iujTH) 2 in the expression for I AAG (| 
evaluated the integral 

yexp(-x - y) 



(C9) 
we exactly 



dxdydz 



(y + z)y/xy + xz + yz 



7T 

2 



(see eq. (5.25d) of Ref. 13) which was not calculated by AAG. 

The high frequency asymptotics of the functions I x , I 2 , J3 , 1 4 and J5 coincide with the results (C9), while the 
asymptotic results for ig derived from the exact expressions (C8) in the limit loth 1 takes the form 

- 2 1 



1 H 



32D (iiOT H ) 2 



We observe a difference in the numerical prefactor in the exact expression and in the term I AAG (CS) found in Ref. 

H 

A much more important problem is, however, not in this numerical discrepancy, but rather in the f act t hat AAG 
evaluated the integrals only in the limit loth 3> 1 and did not study the behavior of the integrals ( |C7| ) at lower 
frequencies loth ^ 1 at all. Note that the low frequency behavior of these integrals is crucially important because it 
determines the dependence of the weak localization correction on th at sufficiently low— temperatures Tth <j 1. [Let 
us remind the reader that here we are discussing only the limit of strong magnetic fieldsE 3 ] th <C t v beyond which any 
perturbative (in the interaction) calculation of the weak localization correction is meaningless. Clearly, in this limit 
the condition Tth ^ 1 is co mpa tible with the condition Tt v 3> 1.] 

It is easy to observe from (C8) that at low frequencies the integrals Ii, I3, I4, I§ behave as cx l/y/iutTH, the integral 
I2 tends to a constant at ui — > 0, while the integral Iq behaves as Iq cx 1/ [iloth) 3, ^ 2 ■ It implies that the contribution of 
the integrals to the conductivity correction (C6) does not diverge as temperature goes to zero but stays finite 

for any finite th even at T — 0. The contribution of the integral Iq diverges as 1/VT for low T. Note, that in the 
high frequency limit (C9) all the integrals show the same asymptotic behavior, and one could naively conclude that 
all these integrals are of the same origin. In reality, however, only (a part of) the integral Iq can be interpreted as the 
effect of interaction on the weak localization, while all other terms actually contribute to dephasing, at least for not 
very high temperatures Tth < 1. 

Let us split the function Iq (C8) into two terms which give respectively divergent and convergent contributions to 
the integral of I§(lo) over to at low frequencies 
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We will now treat only the first (divergent) term in this equation, while the second (convergent) term will be added 
to the integrals Ix,..,Iq. The divergent term yields the following contribution: 
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where 



is the interaction correction to the conductivity, and 

6<twl = VD^ (C13) 

is the weak localization correction. Within the validity range of the perturbation theory in the interaction, the equation 



(Oil) is valid in the whole temperature interval fro m Tt h <C 1 to Tth 2> 1. The physical origin of the correction (CI 1) 
is quite transparent. Indeed, the interaction term ( |C12| ) contains the photon propagator C R oc 1/a. If one replaces the 
total conductivity by the sum of the Drude conductivity and the weak localization correction, one will immediately 
observe that in the first order in SawL the interaction correction will be transformed as Sac ~ ► Saci^ — Saw l I a). 



Thus we arrive at (Cll). 

This simple consideration clarifies the origin of the correction ( |C11 ). It can be interpreted as the effect of weak 



localization on the interaction correction. The same conclusion follows if one considers the diagrams contributing to 
the integral Iq. They are just the diagrams which yield the interaction c orrec tion multiplied by the weak localization 
correction to the photon propagator. In other words, the contribution ( |Cll| ) originates from the second term in the 
right-hand side of eq. (5.22) in Ref. |l3|, while all the remaining contributions come from the first-term, i.e. from the 



unrenormalized photon propagator. Collecting all these contributions we express the final resultli 3 ] for Sa^^ L (C6) 
as a sum of two terms 



?>°cwl = bacwL + Sa-\ (C14) 



where SacwL is defined in (Cll) and 



~ 1 303/2) ^ <IVDTE. Tt H >1. 




2ttT \ 2 



can be derived from the eq. (C6) combined with (|C 

We observe that for Tth 3^ 1 the contribution proportional to th dominates and determines the behavior of 
<5<7 ? ' ? as a function of tjj. This contribution cannot be interpreted as an interaction correction since it grows faster 
than ^Jth with increasing th- For quasi- Id systems all contributions to the conductance growing faster than y/m 
should be interpreted as dephasing terms (see also our discussion in Section 2). Thus - although the correct and 
complete description of the interaction-induced decoherence cannot be obtained from the perturbation theory in the 
interaction - the conclusion about the existence of temperature-independent dephasing in disordered contLuctors at 
low T (determined by the range of frequencies u > T) follows already from the perturbative calculation^ This is 
in contrast with the qualitative arguments about the absence of zero-temperature dephasing presented in the same 
paper. 

It also follows from the above consideration that the suggested by AAG splitting of the weak localization correction 



into two terms (CI) can hardly be justified in both limits Tth 3> 1 an( i Tth <C 1 even within the framework of 
their perturbation calculation^. Indeed, for Tth 3> 1 the total exp ression for the weak localization correction Sa^wL 



derived by AAG and defined here in eqs. (CI), (|C2|), (C6) and (CI) can be expressed as a series expansion in 1/^JTth 



In this limit Sa{jffi L is vanishingly small (it contributes only to the next after the subleading term in (|C16| ) and can 
be safely disregarded. Moreover, for Tth ^S> 1 exactly the same term ~ &< j cw'l ^ \/ t h/T (with a slightly different 
numerical prefactor) is contained in the expression for Saf^ (|C5|). Also due to this fact there are no reasons to 
distinguish the term 5a^y^ L fr om t he remaining contribution to Soyfjf in the above limit. 



For Ttr ~ 1 both terms in (CI) are of the same order, and therefore their separation is not possible. Finally, in 
the limit Tth <C 1 the result for Sa^f^ reads 

Sa*i G = Sa C wL + --^!^(l + 0(VTT H ~)). (C17) 

7T a\ Z7T 
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The origin of the term SacwL oc ypra /T ( |C11 ) was clarified above. It does not describe dephasing and it is purely a 
matter of convention whether to include this term into the perturbative weak localization or interaction corrections. 



The second term in (C17) is represented as a series expansion in powers of ^Tth, it describes dephasing and remains 
finite even at T = 0. We emphasize again that in the limit Tt r *C 1 the leading "dephasing- type" .-contribution to 
the perturbative weak localization correction comes from 5<Jcw*l an d n °t from the "dephasing" terrrO Sa^fa which 



only contributes to the higher order terms (~ (Tth) 4 ) of the expansion of Sa^f^ in powers of Tth (cf. eq. (C5)). 
Terms of the same order are also contained in the expression for ^cr^WL- Thus the splitting ( |Cl| ) is not justified also 
in the limit Tth <C 1. 



APPENDIX D: 



For .-reference purposes we present some rigorous expressions obtained within the exactly solvable Caldeira-Leggett 
modeB for a quantum particle with coordinate x interacting with an infinite bath of oscillators. The time evolution 
of the density matrix of such a particle is defined by eq. (|l|) and the kernel J is given by the following path integral: 

X lf X 2 f 

J(t,x 1 f,x 2 f,x 1 i,x 2l ) = J Vxi(t) J T>x 2 (t)eyip (iS [xi(t)] ~ iS [x 2 (t)} - iS R [x 1 (t),X2(t)] - S I [x 1 (t),x 2 (t)}) , (Dl) 

Xn X 2 i 



where Sq[x] = J dt' (mx 2 /2) is the action of a free particle. The interaction part of the action has the following form 
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t ti 
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ai{h - t 2 )xi(ti)x 2 (t 2 ) - ai(h - t 2 )x 2 (ti)xi(t 2 )} . 



The kernels ai and a R are given by the following integrals: 

a R (t) =-% J — I{w)e 
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cti(t) — / — I(uj) coth— -e 
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where I(co) is the spectral density of the oscillators. This function can be arbitrary, but we will consider here only the 
case of Ohmic dissipation, I(oj) — r)LU0(oj c — |w|), with uj c being the high cutoff frequency. Th is spectrum is the most 
relevant in view of comparison to the disordered metal. The last term in the action S R (D2) compensates potential 
renormalization cau sed by the interaction and maintains the translational invariance of the system. 

The kernel J (Dl) can be found exactly as the integrals over the coordinates are Gaussian. One finds (see e.g 



J(t,Xif,X 2 f,Xu,X 2i ) = 
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where 
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g x {t) = ^JdsJds' G(s - s ) ~Tt + In 
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APPENDIX E: 

As it was already pointed out, the authors of several recent papers0~0 arrived at conclusions different from our 
and, moreover, argued that our approach! is not-eqraect. On top of the standard arguments (which have been already 
discussed in the bulk of this paper) the authorala~ll3 suggested various additional reasons which could invalidate our 
analysis. In view of that we feel it will be appropriate to address the arguments presented in the above papers. We 
believe that this Appendix can be useful for the reader who would like to follow the details of the discussion around 
the problem of quantuni-decoherence at low temperatures. 

Cohen and Imry (CI)ll3 proceeded within the Feynman- Vernon influence functional formalism and found that the 
interference of any pair of time reversed paths is suppressed due to interaction with an effective environment even 
at T =p.Q, This result (eq. (12) of Ref. |l5|) implies-ajjapnzero decoherence rate at T — and is consistent with our 
resultajQ as well with the results of other authorsE£lc3 obtained within the framework of the CL model. However, 
CI argued that the saddle-point approximation they used "cannot be trusted" at low temperatures, and the T- 
independent contribution to the dephasing rate "should be omitted" . The CPs arguments in favor of this conclusion 
follow the line of reasoning according to which a particle with the energy ~ T cannot excite environmental modes 
with energies exceeding T. 

The argumentstJ can not be accepted. Indeed, if CI do not trust the saddle-point approximation which gives 
nonzero dephasing down to T — 0, this could only imply that one should analyze the role of fluctuations around the 
saddle points. The contribution of non-saddle-point paths may only yield further suppression of quantum coherence 
simply because the relevant saddle points provide the (local) minimum for the action. The imaginary part of the 
effective action is positive Si > for all paths except for pairs of exactly equal paths in which case Si = 0. [The 
latter paths do not contribute to the dephasing rate.] Thus the saddle point approximation may only underestimate 
the dephasing rate. 

We can add that in our problem the applicability of the saddle point approximation cannot depend on temperature. 
This is particularly clear in the weak interaction limit. In this case the saddle point paths are determined only by 
the "noninteracting" parts of the action which do not depend on temperature at all. The "interacting" contribution 
can then be treated perturbatively in the exponent (cf. e.g. Sec. 5C), and this is a completely legitimate procedure 
controlled by a small parameter 1/pfI <C 1 in the case of disordered metals. At sufficiently high T this procedure 
yields the AAK's resultH. If one accepts the saddle-point approximation at higher T, one should accept it also at lower 
T: the saddle point paths are the same and the "interacting" contribution to the exponent may only decrease with 
decreasing T. The saddle point approximation which treats the interaction term perturbatively in the exponent (and 
thus yields energy conservation on the saddle point paths) is well known and was frequently used for the Feynman- 



VernpB influence functionals, see e.g. Refs. p4]35 



also mentioned the simple perturbative results which - according to them - yield zero decoherence rate at 
T = 0. As it was already explained above, (i) the problem is essentially nonperturbative and no information about the 
dephasing time t v can be extracted from perturbation theory in the interaction (Sec. 2) and (ii) at low temperatures 
the golden rule approximation yields incorrect results even for perturbative terms (Sec. 4B). Eqs. (14), (15) of Re f. [l5| 
can be obtained from the equation above eq. (14) of the same paper only within the golden rule approximation Q65|). 



Exact calculation leads to an additional term (cf. eqs. (67) and (]l 1 8 ) ) which survives even at T = and diverges at 



large times. Unfortunately this term is missing in eq. (14) of Ref. |15 
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The same term is missing also in eq. (1) of the paper by Imry, Fukuyama and Schwab (IFS)0. Already this fact 
as well as our analysis of Sec. 2 invalidates the "proof that zero-point motion does not dephase" presented by IFS. 
Again, their eqs. (1), (2) are obtained within the same golden rule approximation. Since IFS allow for a general form 
of the interaction matrix elements, it should be legitimate to take ones e.g. for the CL model and to substitute them 
into eq. (1) of Ref. [T| Then one would immediately arrive at the conclusion that no decoherence occurs in the CL 
model at T = 0. Thus the proof of Ref. [l?] explicitly contradicts to the results obtained within the exactly solvable 
Caldeira-Leggett model. 

IFS also argued that dephasing cannot occur at T — > in equilibrium because "in that limit neither the electron 
nor the environment has any energy to exchange". It is well known, however, that the energy of a subsystem (an 
electron in our case) interacting with other subsystem (other electrons) is not conserved due to interaction. The energy 
exchange between different subsystems of a closed system is, of course, possible even at T = 0. In the presence of 
interaction with any other quantum degrees of freedom an electron can be described only by the density matrixE2l and 
is (obviously) not in its nonintera&ting ground state. Therefore it always has energy to exchange. The same is true for 
the environment: it is well knownM that an imaginary part of the dielectric susceptibility for a large-system Ime does 
not vanish even at T = 0, thus implying the possibility of energy exchange. The above argumentO disregards the 
interaction term in the Hamiltonian, and this is again nothing but the golden rule approximation (see also Appendix 
B of Ref. [li] for further discussion of this point). Due to interaction the ground state energy of the total system is 
different from (in our case larger than) the sum of energies of its noninteracting parts, and the energy exchange is 
always possible. PH., 

IFS also pointed out that our resultsM are "in disagreement with experiments" of Ref. |3(| The experimental results 
for the dephasing time t v were reported only in one (the second) out of three papers under Ref. [36|. In this experiment 
it was found that at sufficiently low T (below ~ 1 v 3 K) "the dependence (T) is flattened out". The comparison 
between our results for r v and the experimental findings of the second Ref. |3^ was carried out in Ref. ^. An excellent 
agreement was revealed for all three samples studied in-the second paperE3, see Fig. 3 of Ref. [?]. Therefore, it remains 
unclear which of our results was meant by the authoralU to be in disagreement with experimentsE3. „ 

Since this point was not clarified by IFS we can try to conjecture that they actually interpreted our conclusion^ 
that Coulomb interaction in weakly disordered quasi-ld metallic (many channel) conductors (described within the 
standard Drude model) precludes the Thouless crossover into an insulating state as contradicting to a rapid growth 
of the wire resistance with decreasing temperature detected in Ref. ^ below T ~ 1 K. If so, we can only point out 
that the above conclusion does not contradict to the experimental data but only to their interpretation in terms of 
the Thouless crossover adopted in Ref. The wire conductivity can be represented as a sum of the Drude term o~\ , 
the interaction correction 5o~i n t(T) and the weak localization correction 6o~wl(T, H): 

0-1 + Saint (T) + Sa W L (T, H) . 

Our analysisSi demonstrates that the last of these three terms saturates at low T, and this is in agreement with the 
observations reported in the second Ref. [?6[ However the total resistance may well keep increasing at even lower T 
because of the interaction term do~i n t(T). This scenario indeed precludes the "noninteracting" Thouless crossover, but 
not the crossover into an insulating state due to interaction (e.g. of the Coulomb blockade type). Since in Refs. |,| 
we did not address thp_.term 5o~i n t{T) at all, one can hardly argue about any contradiction between our results and 
the experimental dataEx ■— ■ 

Finally let us turpJrt .the paper by Vavilov and AmbegaokarLLa These authors did not employ the "orthodox" golden 
rule approximationtj'EZfE3 but rather attempted to.-analyze the problem by means of a high temperature expansion. 
They also presented a critical analysis of our paperO (see Appendix C of Ref. |l^) . 

VA questioned the validity of our procedure which amounts to deriving a dephasing time t v only from the terms 
in our expression for the effective actiontj which dominate at sufficiently long times. According to VA the dropped 
terms might be important at times t ~ t v (where Si{t v ) ~ I) and due to that the result for t v could be different 
from ourstl Even without making any calculation one can realize that the contribution of these dropped terms, if 
important, could only make the dephasing time shorter than that found in our papera. Indeed, if one assumes that 
taking all terms into account one would obtain a longer dephasing time r' 3> t v , one would immediately arrive at a 
contradiction with the fact (acknowledged by VA) that the dropped terms are unimportant at times t ^> t v : at least 
at times t ~ t' v (and, hence, at t ^> t v ) these dropped terms should still be significant. Since this is not the case, by 
neglecting the above terms one actually gets an upper. bound for t v . This is in contrast with the VA's claim that t v 
is parametrically longer than that found in our paperQ at low T. n 

The_analysis presented in Sec. 3A of this paper fully confirms our previous resultsu. All terms of our effective 
actiontj were explicitly taken into account in eq. (|23|). Additional terms in the exponent (see eqs. (p8|), (|3(i|)-(34)) 



indeed appear, but (i) they lead to further suppression of quantum coherence and (ii) for all times t > t v they are 
small as compared to the leading order terms which we already considered beforeO. 
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The problem with "unphysical divergences" discussed in the beginning of Appendix C of Ref. R3| does not exist 
in our analysis either. This can be observed already from the footnote 32 before eq. (76) of Ref. 6f the actual low 
frequency cutoff in this equation is at cj ~ l/t, see also eqs. (46), (47) of Ref. |?j. This cutoff is not imposed by hand, 
but rather follows from the fact that the integral over ui in eq. (76) of Ref. |(] and related formulas is just the long time 
approximation of the sum over discrete Fourier frequencies w„ = 2nn/t. [The cutoff at ~ l/ r ¥= was used e.g. by AAKcl 
and in many other papers, we have followed this approximate and physically correct procedure in eq. (76) of Ref. ^ 
just in order to illustrate the similarity of the low frequency contribution to t v found in Refs. || and |^.] The absence 
of any "unphysical divergences" in the exact expression ( |23|) is also completely transparent. It is hard to understand 
why VA believe that our results cannot be compared to those obtained within the Caldeira-Leggett model. A close 
similarity between both problems is obvious from our analysis presented in Sec. 5. 

Another problem of VA wit h our analysis has to do with the factor (1 — 2n(p, r)) which appears in the real part Sr 
of our effective action (see eq. (A18)). VA stated that in Ref. ^we "neglected the time dependence of the momentum". 



This is not true. The electron momentum changes its direction after each scattering event, an this fact is explicitly 
taken into account in Ref. || where the classical electron trajectories in a disordered potential were considered. No 
momentum conservation was imposed in our analysis, rather the electron energy conservation on these trajectories 
was used. This approximation is fully justified, since the saddle point paths are determined by the "noninteracting" 
terms in the effective action while "interacting" terms are treated as a perturbation in the exponent. For these saddle 
point paths the factor (1 — 2n(p, r)) (and not pit)) does not depend on time. It is not quite clear why VA suggest to 
consider odd functions of time nLn(t)). Such functions do not contribute to the exponent at all. It is also not clear 
in which context the observation!!^ that for an odd function of time n(p(t)) the real part of the action Sr is of the 
same order as Sj (taken at T = 0) could be important: there is no way how the term iSr can cancel Si for real and 
nonzero Sr and Sj. The contribution of the trajectories with nonzero Sr to the path integral may only be suppressed 
further due to the presence of the oscillating term exp(— iSr). 



Trying to justify their arguments VA presented several equations (eqs. (C12-C15) of Ref. 16) which are somewhat 
reminiscent to ones obtained e.g. in Sec. 4 of this paper. For instance, the combination "coth+tanh" appears in eq. 
(C15) of Ref. H in a correct form (^33|) rather than in the form (|62|) used by some other authors. Unfortunately VA 



did not evaluate their eq. (C15) but just concluded "we see that the contribution of high frequencies is exponentially 
suppressed" . In Sec. 4B we have demonstrated just the opposite: the high frequency contribution is not exponentially 
suppressed even at T — > and, moreover, it leads to the presence of diverging terms already in the first order in 
the interaction (cf. eq. (|67|)). Exactly the same terms are contained in eq. (C15) of Ref. [if] and, hence, the above 
statementEJ is explicitly incorrect. In Ref. [l6| we also failed to find a proof of even stronger VA's statement that 
one can "observe the same type of cancellation of high frequency contributions to the conductivity in any order of 
perturbation theory" . n 

Having observed that none of our resultsa suffers from the VA's "critique" let us come to the analysis of the 
calculation!!^. In^arder to understand the problems with this calculation we have to start from another paper of 
the same authorsEH. VA acknowledged that the papeitll "is simply wrong, because it treats the phase as a single 
particle and looses the physics associated with the exclusion principle" . We cannot agree with the second part of this 
statement. In our opinion, it was a good idea of Ref. ^ to apply the Caldeira-Leggett model in order to qualitatively 
understand the experimental resultsa. It is not the absence of the exclusion principle in this model, it is the calculation 
of Ref. |3?] which is problematic. In particular, we mention that according to Ref. ^ the real part of the action Sr is 
responsible for the low temperature saturation of t v whereas the contribution of its imaginary part Si is not important 
at T — ► 0. A correct calculation demonstrates the opposite: both in the CL model and in the case of a disordered 
metal the term Si is important, while the term Sr does not contribute to the exponent (and thus to the dephasing 
time) at all. The term Sr vanishes on the saddle-point trajectories both with and without the exclusion principle, 
i.e. the latter does not play any crucial role in the problem in question. 

Revising their earlier paper, VA applied the same method, of calculating the path integral as one employed in Ref. 
|37| . Should the term Sr be kept in the exponent, the resultEO would follow again. In order to avoid that VA expanded 
the exponent to the first order in Sr and in the "quantum" part of Si and made use of the perturbation theory 
developed in Ref. [Hs]. Several comments are in order. 

1. Even if the perturbative calculations^ of the weak localization correction in the high temperature limit were 
correct, the correct dephasing time t v could not be extracted from these calculations. The reason for that 
is exactly the same as for the perturbative calculationE£l: the T-independent linear in time perturbative con- 
tributions coming from the Taylor expansion of the exponent and the pre-exponent exactly cancel each other 
in the first order. As it was explained in. Sec. 2 the problem of finding the dephasing time t v is essentially 
nonperturbative and, hence, the methodlilj should fail already at the point where the exponent was (partially) 
expanded in powers of the interaction. 
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2. The perturbative analysiJll cannot simply be repeated for the situation considered by VA because the initial 
(non-perturbed) propagators depend on the classical part of the fluctuating field. This implies that the Fourier 
transformation in time cannot easily be performed, and the whole calculation should be redone from the very 
beginning. 

3. In order to obtain the "classical" part of the weak localization correction VA used a minimization procedure, the 
meaning of which remains unclear. Indeed, let us rewrite eq. (92) of Ref. [l6| for the weak localization correction 
in the presence of interaction (for the sake of transparency we keep the notations adopted in the above paper) 

Aa wl = Ao$ + Aa deph + Aa' deph , (El) 

As compared to eq. (92) of Ref. [If], in (El) we omitted the last (smallest) term Aa cw i which is unimportant 
for our discussion. According to VA the first term Aa^ comes from the "diagram without electron-electron 



interaction with finite dephasing rate 7" . This should imply that the rate 7 in Ref. 16 is not due to the electron- 



electron interaction but rather due to some other dephasing mechanism. E.g. in the presence of the magnetic 
field one has 7 = 1/th and eq. (93) of Ref. |l6| for Acr^ is equivalent to our eq. (^) or to the first equation 
(4.2b) of Ref. |l3|. The second term Ao d eph is equivalent to the first term in eq. (4.13a) of Ref. [II] (or to the 



first term in eq. ( 106 ) of the present paper)) provided one assumes 7 = 1/th- This term already originates 
from the electron-electron interaction (classical Nyquist noise) , but again the decoherence rate 7 is not related 
to this interaction but rather is due to some other mechanism. We also note that the term Aa de ph written in 

the form of eq. (96) of Ref. Eflmakes sense only if it is small Aadeph <C Aa. 



(o) 
ii ■ 



The next step of the authorala is not easy to understand. They minimize the sum of the first two terms in 
(El) with respect to the dephasing rate 7 and in this way they "self-consistently" derive 7 ~ 1/t aak oc T 2 / 3 . 



Unfortunately VA did not clarify (i) why the result (El) should be minimized with respect to 7 and (ii) why the 



rate 7 depends now on the electron-electron interaction though initially it was assumed that 7 is not related to 



the electron-electron interaction at all. [Actually VA minimized only the two first terms in (El). This implies 



that these terms can be of the same order in contradiction to the inequality Aa de ph <C Aa^\. No explanation 
of the above procedure was given in Ref. [l6| 

4. VA found a negative "quantum" correction Aa' deph , see e.g. eqs. (88-89) of Ref. [l6| From this result VA 
concluded that the dephasing ratea "is overestimated, and the quantum correction suppresses it". In other 
words, according to the authorsli3 quantum noise - at least at not very low T - acts "against" the classical 
noise. No physical reasons supporting this strange result was presented in Ref. [l6]. Fortunately there is no need 
to lppk for such reasons. It is obvious e.g. from our eq. (|30|) that the correction to the classical dephasing 
rateo is positive at any relevant temperature, i.e. quantum noise may only enhance the dephasing effect of 
the classical noise. Below we will demonstrate that no conclusion about both the sign and the value of the 
quantum correction to the classical dephasing rate can be drawn from the expression for the leading order high 
temperature correction Aa' deph . 

5. VA's procedure can easily be tested with the aid of the results derived in the present paper. Since the authorsEl 
assume "that the electron-electron interaction is the only mechanism of decoherence" and consider only the high 
temperature limit, in our formulas of Sec. 2-4 we can put 1/th = and keep only the term 




hd{t) = - — d-Tt*'* (E2) 



in the exponent of (||) expanding this exponent to the first order in fi(t) — fi c i{t). We also expand the pre- 
exponent Ai(t) to the first order in the interaction Ai(t) ~ A { ° ] (t) + A^\t), where A { ° ] (t) = l/(2y/nDt). 
This procedure should exactly correspond to the high temperature expansion of Ref. Il6l Then for the weak 



localization correction (El) (which we now denote as 5o~wl) we will obtain 



e 2 VD f dt f m e 2 y/D f dt ( . , „ A\ 1 '(t)\ , m 

*™ = — j ~q e ~ vt [ m - flcl{t) - aW) ) e ■ (E3) 

v 1 w/ 

Here the first integral determines the "classical" part of the weak localization correction while the second int egra l 
gives the "quantum" correction. Both integrals can easily be evaluated. Combining eqs. (|30|), ( |83| ) and ( |105| ) 
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(the last equation allows to evaluate the term A^\t)/A^\t)) and keeping only the leading corrections from the 
high temperature expansion we obtain 



Mt) - Mt) -m a (m*JE t 0. (E4) 

A { °\t) V2 O-l V 7T 

As before, the T-independent linear in time terms contained in fi(t) ( ^0| ) and in A] 1 ' (t)/A^ (t) exactly cancel 
each other. This cancellation illustrates again why the correct dephasing time at low T as well as the leading 
order high temperature correction to the classical de p hasi ng rate cannot be recovered by means of a high 



temperature expansion for the conductance. From (E2)-(E4) we find 



^ f geVim 173 C(l/2) e 2 VZ> 

We conclude that the procedure developed by VA in Ref. [l6] does not provide correct information about the 
interaction- induced dephasing in disordered conductors. 
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